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19 Abstract 

20 

2 1 Three-dimensional reconstructions of six rocket-triggered lightning channels are derived from 

22 stereo photographs. These reconstructed channels are used to infer the behavior of the current in 

2 3 return strokes above the ground from current waveforms measured at the channel base and 

2 4 electric-field-change waveforms measured at a range of 5.2 km for 24 return strokes in these 

2 5 channels. Streak photographs of 14 of the same strokes are analyzed to determine the rise times, 
2 6 propagation speeds, and amplitudes of relative light intensity for comparison with the electrical 

2 7 inferences. Results include the following: 1) The fine structure of the field-change waveforms 

2 8 that were radiated by these subsequent return strokes can be explained, in large part, by channel 

2 9 geometry. 2) The average 10 - 90% rise time of the stroke current increased by about a factor of 

30 seven in our sample, from an observed 0.31 ± 0.17 us at the surface to an inferred 2.2 ± 0.5 ps at 

31 1 km path length above the surface. 3) The three-dimensional propagation speed of the current 

32 front averaged 1.80 ± 0.24 X 10 8 m/s over channel lengths typically greater than 1 km. 4) 

3 3 Assuming that the measured current was entirely due to the return stroke forced an unreasonably 

3 4 large and abrupt reduction in inferred current amplitude over the first few tens of meters above 

3 5 the surface, especially in cases when the leader was bright relative to its stroke. Therefore, a 

3 6 significant fraction of the current at the surface was probably due to the leader, at least in such 

37 cases. 5) Peak return-stroke currents decreased by approximately 37 ± 12% from 100 m to 1 km 

38 of path length above the surface. Because of uncertainty about how to partition the measured 

3 9 current between leader and return stroke, we are unable to infer the variation of current 

4 0 amplitude near the ground. 



Introduction 


41 

42 

4 3 The return strokes in cloud-to-ground lightning flashes are the most powerful known 

4 4 lightning processes in terms of both local energy dissipation and electromagnetic radiation [e.g . ; 

4 5 Guo and Krider, 1982; Willett et ah, 1990], Although the threat that they pose to objects on the 
4 6 ground has now been reasonably well quantified, the peak electric currents, current rise times, 

4 7 and electromagnetic-field intensities above the surface are not known. These quantities have 
4 8 practical importance to the hardening of aircraft and missiles, which can trigger lightning 

4 9 discharges in flight [Mazur et ah, 1984; Christian et ah, 1989], since the flashes so triggered are 

5 0 likely to contain return strokes if initiated close enough to the ground. "In testing to determine 

5 1 the immunity of a system to a direct lightning strike... the flash current is generally that observed 

52 at the base of a severe flash to ground even though aircraft and space vehicles in flight will likely 

53 encounter the smaller currents associated either with the upper portion of return strokes or with 

54 various in-cloud currents from either ground or cloud discharges" [Uman, 1988]. 

55 In addition to this engineering application, information about the evolution of speed, 

5 6 shape, and amplitude of the upward-propagating current waveform is needed to test, and to guide 

57 the further development of, physical models of the return stroke [e.g., Strawe, 1979; Mattos and 

58 Christopoulos, 1990]. Such models give predictions of these quantities that (except for 

5 9 propagation speed) have yet to be compared to observations above the surface. 

6 0 Direct measurements of the current waveforms in both natural and rocket-triggered return 

61 strokes exist, but only at the channel base [Berger et al, 1975; Eriksson, 1978; Garbagnati and 

62 Lo Pipero, 1982; Leteinturier et al., 1990, 1991; Fisher et ah, 1993; Depasse, 1994, Crawford, 

63 1998; Uman et ah, 2000; Schoene et al., 2003]. Time-resolved photographic and photo-electric 
6 4 measurements that show the evolution of luminosity along the visible channel on microsecond 

65 temporal, and several-meter spatial, scales are also available [Schonland et ah, 1935; Idone and 

66 Orville, 1982; Jordan and Uman, 1983; Idone et ah, 1984; Mach and Rust, 1989; Jordan et ah, 

67 1997; Wang et ah, 1999]. These optical measurements imply that the propagation speed and the 

6 8 peak amplitude of return-stroke current wavefonns both decrease with increasing height above 
69 the surface, while the current rise time increases with height. 

7 0 Unfortunately, the relation between optical emission per unit channel length and 

7 1 instantaneous current is not well known under the conditions found in lightning return strokes. 
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72 Colvin et ah [1987] found an approximately linear relation between instantaneous current and 
7 3 channel brightness in oscillatory laboratory discharges having current periods of about 2 ms. 

7 4 These discharges were intended to simulate "nuclear lightning," however, so their time scales of 

7 5 were slow enough that pressure equilibrium was well established within the current-carrying 

7 6 channels. Murphy et ah [1986, Figure 13] made similar measurements on a laboratory discharge 

7 7 ringing with about 30 ps period and having a current rise rate of approximately 10 kA/ps. (Note 

7 8 that this current derivative is roughly one order of magnitude smaller than those typically 

7 9 measured at the channel base in rocket-triggered return strokes [Leteinturier et 'al. , 1990, 1991; 

8 0 Schoene et ah, 2003] but might be comparable to those occurring in lightning channels a few 

81 hundred meters above the ground.) Murphy et ah [1986] concluded, "This discharge rings on a 

82 timescale too short for pressure equilibrium to be established with the surrounding atmosphere, 

83 thus its brightness is not a simple function of the discharge current." Gomes and Cooray [1998] 

8 4 measured the optical emissions from the entire length of laboratory sparks in 250 mm and 500 

8 5 mm gaps. They found that, for a current rise time of 1 . 1 ps, the optical rise time remained near 

8 6 1 .2 ps over a range of peak currents from 0.6 to 3.5 kA; that optical and current rise times were 

8 7 nearly proportional over a range of 0.3 to 1 5 ps (peak current varied with rise time); and that 

8 8 peak optical emission was approximately linearly related to peak current, although with a large 

8 9 zero intercept that eliminated any proportionality between these two parameters. 

90 In measurements on rocket-triggered lightning, Idone and Orville [1985] reported 

9 1 approximate proportionality between peak current measured at the ground and peak relative 

92 luminosity about 50 m above ground for 20 strokes in one flash and 17 strokes in another. Wang 

93 et al. [2005] found proportionality between instantaneous current and instantaneous relative 

94 intensity in the lowest 3.6 m of channel during the rapidly rising portion only of four triggered 

95 strokes. Nevertheless, the contradictory nature of the laboratory results to date makes it 

9 6 premature to deduce current behavior directly from optical-luminosity measurements. 

97 (Note further that rocket-triggered lightning differs from natural cloud- to-ground 

98 lightning in that the former has no analog to the natural first return stroke [e.g., Uman, 1987, 

9 9 Section 12.5], Le Vine et ah [1989] have argued the similarity between rocket-triggered and 

100 natural subsequent strokes, based on wide-band recordings of their radiation fields, but the 

101 current waveforms in natural first strokes are expected to be different for various reasons [see 
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102 review and discussion by Willett et al., 1995], Thus, the results in this paper may not apply to 

103 natural first return strokes.) 

104 There are numerous, semi-empirical, "engineering" models of return-stroke currents in 

105 the literature [e.g,, Bruce and Golde, 1941; Uman and McLain, 1970; Lin et al., 1980; Hubert, 

10 6 1985; Heidler, 1985; Diendorfer and Uman, 1990], These and other model variations have been 

107 reviewed and extended by Nucci et al. [1990], Thottappillil and Uman [1993], and Rakov and 

108 Uman [1998]. Many of these models have been tested against observed channel-base-current 

109 and remote electromagnetic-field waveforms [Lin et al., 1980; Thottappillil and Uman, 1993; 

110 Rakov and Uman, 1998], but only assuming a straight, vertical lightning channel and a "realistic" 

111 profile of propagation speed. Le Vine and Willett [1995] have presented evidence, however, that 

112 the channel morphology plays an important role in deteimining the structure of the radiation- 

113 field waveform, even for subsequent return strokes. 

114 Master et al. [1981] computed the electromagnetic fields to be expected aloft from a 

115 nearby return stroke, based on a modification of the model of Lin et al. [1980], but they could 

116 find no measurements with which to compare their results. More recently, Reazer et al. [1987], 

117 Mazur et al. [1990], and Mazur and Moreau [1992] have presented in-sitn observations of a few 

118 events in direct lightning strikes to aircraft that they believe to have been return strokes, but the 

119 evidence for this claim is problematic. Reazer et al. [1987] found peak currents of 1 -4 kA and 

120 rise times around 2 ps in these events, in agreement with the general expectation that return 

121 stroke peak currents should be smaller, and current rise times should be longer, aloft than at the 

122 ground. Mazur et al. [1990] and Mazur and Moreau [1992] did not report peak currents or rise 

123 times, but the latter authors state, "...return stroke currents at flight altitudes are much smaller in 

124 amplitude than those measured on the ground, and are usually smaller than current pulses of dart 

125 leaders and recoil streamers. This observation strongly indicates the need for reexamining the 
12 6 threat to aircraft from return strokes. " 

127 

128 

12 9 Data and Approach 

130 

131 The present paper exploits an existing set of channel-base-current recordings, stereo still 

132 photographs, streak photographs, and remote electric-field waveforms of rocket-triggered 
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133 lightning return strokes to deduce quantitative features of their current waveforms above the 

134 ground. This uniquely comprehensive data set was obtained during the summer of 1987 at the 

135 NASA Kennedy Space Center in Florida and has been described in detail by Willett et al. [1989]. 

13 6 A novel feature of these data is the stereo pairs of still photographs, which have enabled 

137 piecewise-linear reconstruction of the actual three-dimensional geometry of six cloud-to-ground 

138 channels, as outlined previously by Willett and Le Vine [1995], (See Appendix A for further 

139 details.) The reconstructed channels are believed accurate to a few tens of meters or better. 

14 0 They are smooth and approximately straight over the lowest few hundred meters, where the 

141 lightning followed the triggering wires, and tortuous above. Both current and electric-field- 

142 change waveforms were recorded during a total of 24 return strokes in these channels, and 

143 stroke-propagation speeds were also measured from the streak photographs for 14 of these 
14 4 strokes in five channels. 

145 It has been shown theoretically [Hill, 1969; Le Vine and Meneghini, 1978; Le Vine and 

14 6 Kao, 1988; Cooray and Orville, 1990; Vecchi et al., 1994] that kinks and bends in a lightning 
147 channel should produce signatures in the return-stroke radiation field. For example, if an 
14 8 unchanging current waveform were to propagate at constant speed — the "transmission-line 

14 9 model" (TLM) of Uman and McLain [1970] - up such a tortuous channel, each successive kink 

150 would radiate a facsimile of the current waveform with an amplitude, polarity, and time delay 

151 determined by the geometry of that kink relative to the observer [ e.g ., Le Vine and Willett, 

152 1992], Thus, it is tempting to conclude that knowledge of the channel geometiy would permit 

153 deconvolution of the radiated field to yield the spatial evolution of the current waveform along 

154 the channel. 

155 Mathematically, this inverse problem cannot be solved, of course, since the return-stroke 

15 6 current is, in general, a function of both time and path length along the lightning channel, 

157 whereas the field waveform observed at a single location is a function of time alone. Even if the 

158 current were constrained to reduce its dimensionality from two to one, the solution might not be 

159 unique. Here the inverse problem is avoided by solving the forward problem with various 

160 assumptions about the current distribution and comparing the results to observation. Preliminary 

161 results of this approach have been presented by Willett et al. [1989] and by Willett and Le Vine 

162 [1996,2002,2003]. 

163 
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Calculation of the Field Change 

Le Vine and Meneghini [1983, Eq. 26] derived the following equation in vector notation 
(which has been repeated as Le Vine and Willett [1992, Eq. 1]) for the total electric field due to a 
TLM current waveform, Irufft - 4] - l[r’ - r a ' ]/v), propagating on an arbitrarily located, 
arbitrarily oriented, short, linear channel segment in free space. Here v is the propagation speed 
of the current waveform, I is a unit vector in the direction of propagation along the channel 
segment (also taken to be the positive direction for current flow), and 4 is the retarded time that 
the onset of the waveform, Itlm( 0), arrives at the origin of the channel segment, r a (This 
"retarded" time is delayed by the interval, \r a ' - r\!c, that is required for information to propagate 
at the speed of light, c, from the origin of the segment to the observer's location, /*.) 


E (r, t) = - 


do c 


do 

4 7T 


J segment 


[Itim] {1- (1 • VR) vr 


ds 


R 


4 7T 


do c 
4 7T 


J segment 


J segment 


[Itim] {1-3 (1 • VR) vr} 


ds 


R 2 


{ [Itim] {l - 3 (l • VR) VR} dt '} 


ds ' 
R 3 


( 1 ) 


The integrals in (1) are line integrals that are evaluated along the channel segment. The 
integrations are to be done in the "primed" coordinate system, where ds' denotes a differential 
length along the segment, and V7? is effectively a unit vector pointing from the source point, r', 
to the observer along the separation distance, R = \r- r'\. The square brackets in the integrands in 
Equation 1 only, ([function]), denote the retarded value of the argument of the enclosed function, 
and the dot above the brackets denotes a derivative with respect to the single argument of that 
function. 

Starting from the first term in (1), Le Vine and Willett [1992, Eq. B3] derived a simple 
formula for the electric radiation (i. e. , "far") field, measured at the surface of an infinite, 
horizontal, conducting plane, due to a current-carrying channel segment above that plane: 
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(r, t) = -MZ- 

2 7T Rq 


sin 


( 2 ) 
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212 

213 

214 


/A A . ^ 

(©) e • l z 

— [I tim (t - t a ) - I TIM (t ~ tb) ] 

[l-^ {1 - VRo) ] 


In this formula we have assumed a spherical coordinate system, (R 0 , 6, <p), that is centered on the 
channel segment but has its symmetry axis, z , oriented vertically — perpendicular to the 
conducting plane. (Thus the zenith angle, 6, is greater than 90° for sources above the plane.) R 0 
is the distance from the center of the segment aloft to the observer on the ground plane, fio is the 
magnetic permeability of free space, tb is the retarded time that the onset of the current 
waveform arrives at the termination of the channel segment, r h ', and I tla-^ argument) is the 
functional form of the current at the center of that segment. 

The "induction" and "static" components of the complete electric field due to the same 
channel segment can be readily computed from the second and third terms of (1), respectively. 
They and are given in the same notation as (2) by Equations 3 and 4: 


Eind (r, t) 


do C r b ' 


2 7 X Ro 2 
• Itlm (t 


-1 — [sin (0) (© • l) + 2 cos (0) (l ■ VRo) ] 
t a + tb \ 

9 I 


( 3 ) 


Estat ( r r t ) — 


do C r b ' -r a ' 2 

2 71 Ro 3 

t , 

Itlm (t ’ - 


[sin (0) [© • l) +2 cos (0) (l • VRo) ] 
dt ■ 

9 / 


( 4 ) 


(We have used the average retarded arrival time on the segment in question in order to minimize 
systematic errors in Equations 3 and 4.) 

As with (1), Equations 2-4 require the current to obey the TLM, but only within each 
segment. Thus, any desired current distribution can be approximated by dividing the entire 
light nin g channel into short, linear segments and varying the current wavefoim and propagation 
speed appropriately between successive segments. These three equations are used here to 
calculate the total field change that would be produced at the location of the electric-field-change 
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sensor by various assumed current distributions over all of the segments that make up each 
piecewise-linear, reconstructed channel. Further details on this procedure are given in Appendix 
B. 


Formal Assumptions about the R/S Current Distribution 

As indicated above, we expect both the propagation speed and the peak amplitude of 
return-stroke current waveforms to decrease with increasing height above the surface, while the 
current rise time should increase with height. Thus we adopt a "generalized TLM" current 
distribution similar to that of Cooray and Orville [1990], The current as a function of time and 
position on the lightning channel is given by 


i (t, s) = a (s) I (ti, s) 


(5) 


where t is time measured from stroke onset at the surface, 5 is path length measured upward from 
the surface along the tortuous channel, and a(s) is an amplitude factor that allows for permanent 
charge deposition along the channel. I(tj,s) = 0 for tj < 0, and the current onset is assumed to 
propagate monotonically upward with position-dependent "TLM velocity," v(s). t } accounts for 
the resulting propagation delay as a function of s : 


The current waveform that was measured at the surface, i 0 (t), is smoothed by 
convolution, 


rti 

I (ti, s) = i 0 (ti - t ' ) K (t ' , s) dt ' , ti > 0 (7) 

Jo 

where the limits of integration result from the requirement that both i 0 (t) and K(t, s) vanish for t < 
0. We use the following form for our "causal" smoothing kernel: 
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K (t \ s) 


■c (S) t 1 2 

2 r 3 (s) 


> 0 


( 8 ) 


which has been normalized to conserve the total charge that is transported up the channel by any 
given current pulse. Therefore, although the peak amplitude of I(ti,s) typically decreases with 
increasing s as a result of increased smoothing, a(s) scales the total charge passing any 5. The 
"equivalent width" (the width of a rectangle with the same peak magnitude and total area) of this 
convolution kernel is [r(s) e 2 ]/2. The forms chosen for the various parameters defined above are 
as follows: 


a ( S ) - [ a m i n i + ( 1 - ajninl ) e S,/ al ] [ a min2 + ( 1 ~ a rain2 ) e ^ a2 1 ( 9 ) 


V (s) = v min + (v max - Vmin) e s/I Vl + v^pi e s/L v2 


( 10 ) 


r (s) = r max (1 -e“ s/L t) 


( 11 ) 


where ci m \ n u Lah La 2 ? Vmin> Vmaxi Vamph L v i, L v 2 , T, na x, and L t are constants. 

Note that the generalized-TLM form that is adopted here (as opposed to the piecewise 
TLM requirement in the previous section) restricts somewhat the possible current distributions 
along the return-stroke channel. Nevertheless, the present assumptions allow us to explore the 
current variations that we expect, based on optical observations, while introducing a manageable 
number of free parameters. An understanding of the effects of these various parameters can be 
obtained by examining the examples of I(ti,s) in Figure 1 and of v(s) and a(s) in Figures 5 and 7, 
respectively, below. (The corresponding free-parameter values are given in the figure captions.) 

Eleven free parameters may seem like a lot to fit a given observed electric-field-change 
waveform. Note, however, that we are really just allowing the current amplitude to vary on two 
different height scales (often required to fit the peak radiation field, as illustrated below), the 
TLM velocity to vary on two scales as well, and the current rise time to vary on one height scale. 
In a later section, a physical explanation is offered for on e of the amplitude scales, probably 
rendering it moot. The second TLM- velocity scale is required for only a few strokes. Finally, it 
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is not the values of the individual free parameters themselves, but rather the characters of the 
27 6 variations of current amplitude, rise time, and propagation speed with height, that are the real 

277 objectives of this investigation. Therefore, our fitting procedure does not turn out to have as 

27 8 many degrees of freedom as it might seem. 

279 

280 

281 Detailed Example of Fitting Procedure -- Stroke 8732/2 

282 

283 — The approach that is used in this paper to deduce the evolution of return-stroke current — 

2 84 with height comprises the following steps: 1) Guess the current parameters in Equations 9-11. 
285 2) Compute the magnitude of the resulting field change, Eft), at the observing site from 

2 8 6 Equations 2-4, according to Appendix B. 3) Compare Eft) with the measured field change for 

2 87 the same stroke, E m (t). 4) Iterate until a good fit is obtained. We shall see that this procedure 

2 88 gives reasonably definite results, in spite of the apparent ill-posed-ness of the general 
2 89 mathematical inversion problem. As an example, we examine in detail the fitting of stroke 2 in 
2 90 flash 8732 (stroke and flash identifiers as in Willett et al. [1989]). 

291 

2 92 Smoothing and Extrapolation of the Current Waveform 

293 

2 94 The channel-base current for these events was recorded for either 20 ps (flashes 8715 and 

2 95 8717) or 5 ps (flashes 8725 - 8732), including pre-trigger delay. In order to compute Eft) for the 

2 96 entire time interval during which the current onset propagates from bottom to top of the visible 

2 97 channel — typically about 10 ps — it was therefore necessary to extrapolate the measured current 

2 98 waveforms for more than half of the events in our dataset. Furthermore, the smaller-amplitude 

2 9 9 waveforms were rather noisy, leading to (1) uncertainty in the precise time of current onset and 

300 (2) spurious noise in our piecewise-linear calculation of the radiation field via Equation 2, so 

301 some smoothing was beneficial. (If the current waveform is significantly different on adjoining 

302 channel segments, errors are produced in the radiated field. These errors can be minimized 

303 either by smoothing the channel-base current or by using inconveniently short channel segments 

304 near the ground.) 
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305 In practice, the onset of the measured current waveform was truncated at some small 

30 6 magnitude (less than the digitization interval), and then the remaining record was smoothed with 

307 a tapered, 1 1 -point (60 ns FWHM) moving average. This resulted in a relatively smooth ift) 

308 having a definite onset ti m e (t = 0) without appreciably changing its rise time or wave shape. 

309 Next the decaying portion of the recorded waveform, starting well after the peak, was fitted with 

310 the sum of a constant, a linear slope, and (in many of the cases) an exponential decay. The 

311 resulting analytic shape was used to extrapolate ioit) to later times, as needed. This quasi- 

312 objective extrapolation procedure is not regarded as part of the matching of Eft) to E,ft), per se. 

313 In most cases the results are reasonable (see Figures 2 for two examples), but they appeared to 

314 produce artifacts in a few cases, as we shall see later. 

315 

316 Determination of Current Rise Time 

317 

318 The most obvious result of our analysis, first reported by Willett et al. [1989] and by 

319 Willett and Le Vine [1996], is that the current rise time must increase rapidly with height (or 

320 with path length, s) in order for Eft) to resemble E m (t). Thus, the first step in our fitting 

321 procedure was always to adjust r max and L T in Equation 1 1 in order to obtain an Eft) with 

322 approximately the right amount of "fine structure." 

323 Figure 3 illustrates the excessive fine structure that is obtained for stroke 8732/2 if io(f) 

324 (see Figures 2) is allowed the propagate up the reconstructed channel at a constant speed of 1 .7 1 

325 X 1 0 8 m/s - effectively the pure TLM. (Note that this propagation speed along our 3-D channel 
32 6 was chosen to be somewhat faster than the measured 2-D propagation speed of 1.6X10 m/s ± 

327 20% in order to line up the major waveform features in time.) The fine structure on Eft) is 

328 reduced to a reasonable level, as shown by the green curve in Figure 4, when x max = 4.36 ps and 

32 9 L x ~ 7084 m (that is, x(s) is nearly linear with an initial slope of 0.615 ps/km). Unfortunately, the 

330 remaining fine structure is now delayed significantly with respect to that of the observed 

331 waveform, forcing us to adjust the propagation speed. The reason for this delay is that increased 

332 smoothing (via Equations 7 and 8) delays the fast-rising portion of the current waveform by an 

333 increased interval relative to current onset [which, by definition, propagates according to v(,v)] - 

334 see also Figure 1. 

335 
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33 6 Determination of Propagation Speed 

337 

338 The horizontal, solid-green line in Figure 5 illustrates the constant TLM-velocity profile 

33 9 that corresponds to the green curve in Figure 4. In the numerical code it is possible to calculate, 

34 0 at each time step, the propagation speed of the half- amplitude point on the rising portion of the 

341 current front, while accounting for the fact that the current waveform is becoming increasingly 

342 smooth (the rise time is becoming longer) with increasing height. [Later, this calculation also 

n 

343 accounts for a decreasing amplitude factor, a(s), with height.] For v(s) = 1.71 X 10 m/s, the 

344 resulting "front velocity, v eJ fs),” as we call this calculated effective speed, is shown by the green 

345 dots in Figure 5. Notice that, not only is v e fs) significantly lower than v(.s’), but it also increases 
3 4 6 somewhat with height. 

347 In general, we adjusted the TLM velocity profile so that the fine structure of Eft) 

34 8 coincided in time with that of Eft), while striving to keep the front velocity constant or 

34 9 decreasing with height. (The latter was not always possible, however, as shown below.) For 

350 stroke 8732/2 a nearly constant v e fs) was obtained with v min = 2.1 X 10 m/s, v max = 2.6 X 10 

351 m/s, v amp i = 0, and L vl = 500 m, as shown by the red profiles of Figure 5. Here the average front 

352 velocity is almost exactly 1.71 X 10 8 m/s, and the fine structure of the resulting Eft) waveform 

353 is in good temporal agreement with that of Eft), as can be seen from the red curve in Figure 4. 

354 

355 Determination of Amplitude Factor 

356 

357 It turns out to be generally true (at least within the context of Equations 5-11) that the 

358 amplitude of Eft) is directly proportional, not to v(s), but to v e fs). Thus, it is not surprising to 

35 9 find the peak amplitude of the red curve in Figure 4 (corresponding to v e fs) -1.7X10 m/s) to 

360 be about 29% larger than that of the green curve (v e fs) ~ 1 .3 X 10 s m/s). Unfortunately, the red 

361 curve also peaks about 37% higher than the black curve, Eft), indicating that some adjustment 
3 62 of the current amplitude is required. 

363 It is not satisfactoiy in general to eliminate this "over-prediction" of peak field by 

364 adjusting the velocity profile. First, this would normally result in the front velocity's increasing 
3 65 with height, which is considered un-physical. More importantly, it would almost always delay 
366 features of the predicted fine structure relative to those observed. Therefore, the only practical 
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3 67 way to adjust the amplitude of Eft) is to adjust the amplitude factor in Equation 5. In the case of 

368 stroke 8732/2, the parameters, a min i = 0.70, a mi „2 = 0.60, L a] - 60 m, and L a2 - 400 m, result in 

369 good agreement, as shown in Figure 6. 

37 0 The corresponding profile of a(s) is given by the red curve in Figure 7. Notice that two 

371 height scales are usually required — a short one to bring the peak field into agreement with 

372 observation and a relatively long one to tailor the "tail" of the field-change waveform. As 

37 3 mentioned above, non-uniform a(s ) results in a change in the linear charge density on the 

374 channel from before to after the return stroke's passage. This "deposited" charge density is 

375 proportional to -da/ds, as illustrated in the green curve of Figure 7. Notice that an appreciable 

37 6 fraction of the stroke charge is deposited quite close to the surface in this case. This turns out to 

377 be true generally in our dataset and is explained further in a later section. 

378 

379 

380 Fitting Variations and Problem Cases 

381 

382 Adjustment of Eft) Peak Shape Using Velocity Profile 

383 

384 Occasionally it was not possible to match the shape and/or amplitude of the Eft) peak by 
3 85 adjusting a{s) alone. In such cases the initial front velocity could often be increased, at the 

38 6 expense of introducing a second height scale for v(sj, to yield a larger and/or sharper Eft) peak. 

387 Figures 8 and 9 show the best example of this type of velocity adjustment and its results for 

388 stroke 8732/1. 

389 

3 90 Consequences of the Limit, v(s) < c 

391 

392 There is a physical limit on the TLM velocity - the current onset cannot propagate up the 

393 channel faster than the speed of light. As indicated above, we tried to prevent v ej fs) from 

3 94 increasing with height. This was not always possible, however, without violating the constraint 
395 on v(s), especially when a rapid increase in current rise time with height and/or a relatively large 
3 96 v e jfs) was required. The most dramatic example of such behavior is stroke 8728/ 1 1, illustrated 
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3 97 in Figure 10. Nevertheless, we were able to fit Eft) to Eft) quite satisfactorily in this case, as 

3 98 can be seen in Figure 15 a, below. 

399 

4 00 Extremely Fast Decay of Current Amplitude 

401 

4 02 Several cases require both a min i significantly less than unity and L aI < 1 0 m in order to 

4 03 obviate a sharp initial peak on Eft). This implies that a significant fraction of the return-stroke 

4 04 charge is deposited in the lowest ~10 m above the ground in these cases (although the need for a 

4 05 very rapid decrease in a(s) will be examined further below). Stroke 8726/3, for which a minl = 

4 0 6 0.69 and the preferred value of L a i = 5 m, is one of the most extreme cases in this regard. Figure 

4 07 11 compares Ef t) for three different values of L a i with Eft). Evidently the best agreement is for 

4 08 L a i = 5 m. Figure 12 shows the resulting profiles of a(s) and linear charge density on an 

4 0 9 expanded height scale. A possible physical explanation of this and other similar results is 

410 offered in the Discussion section below. 

411 

412 Probable Failure of Current Extrapolation 

413 

414 In a few cases it appears impossible to get a good fit to the latter part of Eft). We have 

415 already seen an example in stroke 8732/1 (Figure 9), where the amplitude of Eft) becomes too 

416 small after the first 6 ps. Since the current for this stroke (see Figure 1 6a) was extrapolated from 

417 a linear fit to less than 3 ps of fairly steeply falling data following the peak, it seems probable 

418 that this phenomenon is caused by poor extrapolation. A more extreme example of this behavior 

419 is stoke 8725/4. The green waveform in Figure 13 shows the original current extrapolation for 

420 this case — simply our standard constant-plus-linear-plus-exponential fit to some 3.8 ps of 

421 descending current record - which falls to zero about 12 ps after onset. The resulting modeled 

422 field change is shown in green in Figure 14 and is seen to disagree with observation after the first 

423 5 ps or so. If we arbitrarily force the current extrapolation to remain high, however (red curve in 

424 Figure 13), the resulting Eft) — red curve in Figure 14, all fitting parameters remaining the same 

425 — agrees much better with Eft). [This is the only case in which we have "tinkered" with the 

426 current extrapolation in order to obtain a better fit to the observed field change.] 
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42 7 [Neglecting the high-frequency noise from the current waveform, notice that there is still 

42 8 too much fine structure in E c (f), relative to E m (f), during the first 5 [is of stroke 8725/4. We 

42 9 could not further increase the current rise time during the early part of this stroke while 

430 maintaining the observed peak field change without violating the constraint that v(s) < c. Such a 

431 rapid increase of x{s) near channel base would significantly reduce v e /j(s) there and, consequently, 

432 limit E c (t).] 

433 

434 

4 35 Summary of Results for All 24 Strokes 

436 

437 Table 1 gives values of the 1 1 free parameters in Equations 9 - 1 1 for each stroke in our 

438 dataset. A few explanatory remarks are in order here. The "Geometry" values in columns 3 and 

439 4 indicate the path length from the surface to the tip of the triggering wire and to the top of the 

4 4 0 reconstructed channel, respectively. "Initial Slope" (column 6) is the only parameter quoted for 

4 41 the kernel time scale, t(s), in cases where L x is much longer than the reconstructed channel, so 

4 42 that the profile of t(s) is effectively linear. Values are omitted for unused parameters in v(s) and 
4 4 3 a(s). Instead of the exponential form specified in Equation 10, a linear decrease of v(s) from an 

4 4 4 initial value of c appeared to be more appropriate for stroke 8726/1, as indicated in column 12. 

4 45 The last col umn of Table 1 gives the retarded time at which the onset of the current waveform 

44 6 arrives at the top of the reconstructed channel ~ the latest observer time, t, at which E c (t) can be 

4 47 computed. 

4 48 Table 2 gives the corresponding values of 10 - 90% current rise time, effective front 

4 4 9 velocity, v e /s), and amplitude factor, a(s), at representative path lengths along the channel in 

450 each case. It is these values that should be considered the results of our matching of E c (t) to 

451 E m (t), whereas the parameters in Table 1 are dependent on the specific form that is assumed for 

452 the current (Equations 7-11). (Recall that the amplitude factor in Equation 5 determines the 

453 total charge passing a given point on the channel.) Note that the average of v eJ ^s) over the entire 
4 54 current-propagation interval — maximum path length reached by the half-amplitude point on the 
4 55 current waveform divided by the actual (not retarded) arrival time of current onset at the channel 

45 6 top — is given in column 8, and that values of v^(l()00 m) are omitted from column 12 for the 

457 few strokes in which this half-amplitude point never reached s = 1000 m. 
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458 Figures 1 5 compare E c (l) and E m (t) for all 24 strokes. The pairs of strokes from three 

459 different flashes in Figure 15a make a familiar point: The field-change waveform s that are 

4 60 produced by strokes within the same flash tend to be more similar than those from different 

4 61 flashes [ e.g Le Vine and Willett, 1995]. This is not always the case, however, as illustrated by 

4 62 Figure 15b, in which a single flash is seen to produce two qualitatively different classes of 

4 63 waveforms. Another example, in which two sharply peaked, relatively simple waveforms 
4 64 coexist with three having more rounded peaks and more pronounced fine structure, is shown in 

4 65 Figure 15c. The remaining eight strokes, all from the same flash, are shown in Figure 1 5d. 

4 66 Notice that we have been able to obtain good agreement in both fine structure and absolute 

4 67 amplitude between modeled and observed field changes over a wide variety of wave shapes. The 
4 68 corresponding extrapolated waveforms of channel-base current 
4 69 are given in Figures 16. 

470 

471 

472 Relative Light Intensity from Streak Photographs 

473 

474 In addition to providing direct measurements of vertically averaged, two-dimensional, 

475 retum-stroke-propagation speed [see Willett et ah, 1989], the streak photographs of 14 of the 24 

4 7 6 strokes in five of our six reconstructed lightning channels were processed to determine relative 

4 77 light intensity (RLI) as a function of time and position. A typical leader/retum-stroke image 

47 8 from our data set is shown in Figure 17. Such images were digitized, calibrated, and converted 
47 9 to time series of RLI at selected vertical levels. The resulting time series were then smoothed 

480 and analyzed to determine both peak amplitude and 10 - 90% rise time of the return strokes, for 

4 81 comparison with similar parameters deduced above for the corresponding current waveforms 

4 82 (see Table 2). This optical analysis is described in detail in Appendix C, and the results are 
483 presented in Table 3. Although we don't attempt to infer current amplitudes or rise times from 
4 84 the optical parameters tabulated here, for the reasons explained in the Introduction, such a 
4 85 comparison is reasonable in view of the approximately linear relation that we found between RLI 

4 8 6 amplitude and peak current over our 14 strokes (like that reported by Idone and Orville [1985]) 

487 and the similarly linear relation between instantaneous RLI and current during the fast-rising 

488 portions of four triggered strokes that was reported by Wang et al. [2005]. 
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4 89 The reconstructed height above ground and the path length along the lightning channel 

4 90 are shown in Table 3 for each measurement level on each streak photograph. (For comparison 

4 91 with the model results in Table 2, the path length must be used.) As detailed in Appendix C, the 

4 92 stroke RLI amplitude is the difference between the upper and lower limits of the fast-rising 
4 93 portion of the corresponding time series. At levels above the first, the lower limit, or "baseline" 

4 94 (also tabulated), can be viewed as an estimate of the residual brightness of the leader channel, 

4 95 just before stroke onset, since the zero of RLI was determined from the "background" film 

4 96 density prior to leader onset. The uncertainty in RLI amplitude is estimated at about ±0.1, based 
4 97 on the typical noise level illustrated in Figure Cl. The uncertainty in rise time is more difficult 

4 98 to estimate, depending as it does on the clarity with which the fast-rising portion is manifest in 

4 99 the RLI waveforms. The error estimates in the last column of the table are based on the scatter 

500 among multiple estimates of rise time, where such existed, as described with respect to Figures 

501 C2. The worst of these uncertainties are about ±23%, although most are much smaller. 

502 The optical parameters in Table 3 have been plotted against the corresponding electrical 

503 estimates in Figures 1 8a and 1 8b. Recall that the model peak current, i p (s), is not identical to the 

504 amplitude factor, a(s), given in Tables 1 and 2. i p (s) is determined by both a(s ) (see Equations 5 

505 and 9) and the convolution smoothing (see Equations 7, 8, and 11), which increases with path 

50 6 length as specified in Table 1 . Therefore, i p (s) has been computed and given in Table 4. (The 

507 values for zero path length are, of course, those directly measured at the surface.) 

508 

509 

510 Discussion 

511 

512 The entire procedure of adjusting the current-model parameters in order to match the 

513 computed E c (t) waveform to the measured E m (t) has been illustrated in previous sections. Not 

514 only is this procedure somewhat subjective, but the various parameters interact with one another 

515 to a greater or lesser extent. We have already seen how the introduction of a current rise time 

516 that increases with height, through x(s) in Equation 8, causes the effective front-propagation 

517 speed, v eJ j(s), to decrease relative to v(s), although this particular interaction does not cause any 

518 ambiguities beyond an occasional encounter with the physical limitation, v(.v) < c. It has also 

519 been pointed out that the amplitude of E c is directly proportional to v e /f, so that speed changes 


18 



520 affect the inferred amplitude factor, a(s). Again, this causes no serious difficulties, since both 

521 the rise time and the effective speed of the current waveform are strongly tied to observed 

522 features of E m (t) . 

523 The amplitude factor, in turn, affects not only the overall amplitude of E c but also the 

524 amplitude of its fine structure. (Compare the red curves in Figures 4 and 6.) Nevertheless, one 

525 cannot normally trade off (for example) increased a(s) against increased r(s), maintaining the 

526 same level of fine structure, because this would make the overall amplitude of E c too large. A 

527 significant concern arises, however, when the extrapolation of the measured current waveform to 

528 later times is uncertain, as illustrated in Figure 13 for stroke 8725/4. In such cases it would be 
52 9 possible to trade off (for example) a higher current extrapolation against decreased a(s) — 

530 maintaining the same overall amplitude of E c but reducing the fine structure — and decreased r(s) 

531 — boosting the fine structure back to the correct level but yielding different inferred current rise 

532 time vs. height. 

533 This last ambiguity exists to some extent in all 14 strokes having short current records 

534 (flashes 8725, 8726, 8728, and 8732), although examination of the measured current and/or 

535 computed field-change waveforms in Figures 1 5 and 1 6 suggests that it might be a significant 

5 3 6 problem in only 8 (strokes 8725/2, 8725/3, 8725/4, 8725/5, 8726/1, 8726/4, 8732/1, and 8732/2). 

537 Since we have no way of assessing the accuracy of the extrapolated currents at late times in these 

538 cases, we cannot accurately evaluate the true uncertainties in the inferred rise time or peak 

539 amplitude of the current at those late times. (This is also why we have refrained from tinkering 

540 with the measured currents in all cases except 8725/4.) The potential uncertainties in these 

541 parameters, primarily at the higher altitudes, due to the extrapolation of measured current 

542 waveforms are therefore ignored in the following discussion. 

543 

544 Rise Time of Current 

545 

54 6 The rapid increase in current rise time with increasing path length, illustrated in Figure 1 

547 for example, is the most robust and inescapable conclusion of the present analysis. We could 

54 8 find no alternative to spreading out the fast-rising current front, at least within the framework of 

54 9 our generalized transmission-line model, to obviate the excessive high-frequency radiation that 
550 would otherwise result from the observed channel tortuosity. The first few rows of Table 5 show 
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551 that the inferred 10 - 90% current rise time, averaged over all 24 strokes in the dataset, increases 

552 by a factor of seven, from 0.3 1 ± 0.17 ps at the surface to 2.2 ± 0.5 ps at 1000 m. (It is worth 

553 mentioning that the results do not change much, especially at the higher altitudes, if we average 

554 only over the 10 strokes with long current records. In fact, the increase from the surface to 1000 

555 m for this subset is nearly a factor of 12, apparently as a result of shorter rise times at the surface 

55 6 for these strokes!) Although there is considerable relative variation from stroke to stroke around 

557 these means, the reduction in coefficient of variation from 55% at the surface to only 22% at 

558 1000 m suggests that much of this is due to differences among the measured current waveforms 

559 at the surface. 

5 60 The accuracy of the individual rise times in Table 2 is difficult to assess, given the 

5 61 qualitative nature of the fitting procedure, but the inferred values appear good to +20% or better. 

5 62 This uncertainty was estimated by varying the rise time of the assumed current and observing the 

5 63 effects on the computed field change in comparison to the observed waveform. In any case, the 

5 64 uncertainty in deduced current rise time is clearly far too small to negate the overall conclusion 

5 65 here. 

5 66 Referring to the right sides of Figures 18, we see that the optical rise times deduced from 

567 the streak photographs are nearly always larger than the inferred current rise times, although the 
5 68 overall trends of these two parameters with increasing path length are roughly parallel in most 

569 cases (notable exceptions being 8725/5 and 8726/3). These results are also summarized in Table 

57 0 5, where the estimated optical rise times have been averaged over all 14 strokes and over all path 

571 lengths within +100 m of the nominal 30, 300, and 1000 m levels. (The mean and standard 

572 deviation of the path lengths included in each of these averages are also given in the table.) 

573 As discussed in the Introduction, there are no direct comparisons in the literature for the 

574 current rise times inferred by our modeling. We can, however, compare our optical rise times 

575 with other optical measurements. Most relevant is the streak-recording analysis of Jordan et al. 

57 6 [1997] on dart leader and return stroke luminosity variation with height. Indeed, the retum- 

577 stroke rise-time variation with height cited by Jordan et al. (1 .5 to 4.0 ps between ground and 1 .4 

57 8 km aloft) is very similar to that presented here in Table 5. Also, Olsen et al. [2004], using a 

57 9 photoelectric technique applied to a triggered flash, provided in their Figure 3 a series of 

580 luminosity traces that have very similar characteristics to ours. Mach and Rust [1989] also 

581 presented some luminosity traces from a photoelectric sensor that are apparently quite consistent 
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5 82 with our derived traces. On the other hand, Wang et al. [1999] reported that the 10 - 90% rise 

5 83 time of RLI in two triggered return strokes increased by factors of 2 and 3.5 over the lowest 40 - 

5 84 45 m of channel, a region rarely imaged in time-resolved photography. 

5 85 Comparing estimated optical and inferred current rise times in Table 5, we see that the 

586 former average about 1 ps longer than the latter at 30 m and 300 m but almost 2 ps longer at 1 

587 km. These differences cannot be fully explained by the approximately 0.5 ps time resolution of 

5 88 the streak camera that was used for these measurements. In addition to the uncertain relationship 

5 89 between luminosity and current in return strokes (see literature review in the Introduction), 

5 90 however, there are two possible reasons for the greater discrepancy at the highest altitudes (latest 

5 91 times) in the present dataset: 1) The uncertainty in our estimates of RLI rise time tends to be 

592 relative (see the last column of Table 3), so that the absolute uncertainty is greater for the larger 

5 93 values. This might be due, in part, to our use of stronger temporal filtering for the slower-rising 

5 94 waveforms, which also tend to have smaller amplitude (lower signal-to-noise). 2) The uncertain 

5 95 current extrapolations in some cases (see the beginning of the Discussion section) result in 

5 96 greater uncertainty in the inferred current rise times (and amplitudes) at late times (high 

5 97 altitudes). 

598 

5 99 Propagation Speed 

600 

601 Table 5 also gives the mean, inferred, effective, current-front-propagation speeds, both at 

6 02 four representative path lengths above the surface and averaged over the full extent of the 

603 reconstructed channel for which this parameter could be computed. The individual inferred 

604 speeds (see Table 2) are fairly consistent from stroke to stroke — they have low a coefficient of 

605 variation around their mean at each level — the only obvious outlier being 8725/4, with an 

60 6 averaged speed of only 0.9 X 10 8 m/s. The means are also fairly uniform with height at about 

607 1.8X10 8 m/s, in spite of the fact that there is a substantial (and apparently unavoidable) increase 

608 in Vefffs) with height in two strokes (8728/10 and 8728/1 1) and an appreciable decrease with 
60 9 height in one other (8726/1). These trends in propagation speed should not be taken very 

610 seriously, as mentioned earlier, but the averaged speeds are probably good to better than ±1 0% — 

611 it's fairly easy to detect and eliminate timing differences between the computed and measured 

612 fine structure. The variation of averaged v e ^-from stroke to stroke is probably also real. 
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613 The mean, two-dimensional, optical propagation speed, V 2 D, over 14 of these same 

614 strokes in five of the six channels is given in Table 5 for comparison, from original data 

615 (averaged over the lowest several hundred meters of channel) reported by Willett et al. [1989]. 

616 Note that our mean, averaged, three-dimensional v^-is about 23% larger than the mean, 

617 averaged, 2-D, optical speed. This difference is not too surprising in the light of previous work 

618 on 2-D vs. 3-D channel lengths by Idone et al. [1984], It can also be seen from Table 3 that the 

619 ratio of the total path length to the vertical height from the surface to optical level nine averages 

62 0 1 .45 + 0.30 over our six reconstructed channels. 

621 Once again there are no direct comparisons for our inferred current-front speeds, but the 

622 originally reported optical two-dimensional mean speed of 1.5 X 10 8 m/s is quite consistent with 

62 3 the results of Idone et al. [1984] and Mach and Rust [1989], both of which had sizeable numbers 

624 of strokes analyzed from triggered flashes over comparable heights. 

625 

62 6 Amplitude of Current 

627 

62 8 The mean, inferred, charge-amplitude multiplier, a(s), can be seen in Table 5 to decrease 

62 9 from unity at the surface (by definition) to 0.68 + 0.08 at 100 m and to 0.49 ± 0.09 at 1000 m. (It 

630 is worth mentioning that the result at the 1 000 m level does change appreciably when we average 

631 only over the 1 0 strokes with long current records. This suggests that a(s) in the upper parts of 

632 the channels might be affected significantly by our extrapolation of the corresponding surface 

633 current records, as mentioned in the sub-section, Rise Time of Current.) The relative variability 

634 of the values for individual strokes at the various levels is quite small, although 8725/3 and 

635 8726/1 might be considered outliers at the lower levels (see Table 2). As before, the uncertainty 

636 in our inferred values of a(s) (and in the peak currents in Table 4, discussed below) was 

637 estimated to be ±10% or better by varying a(s) and observing its effects on the computed field- 

638 change amplitude in comparison to the observed waveform. This low apparent uncertainty 

639 should not be too surprising, since amplitude differences between computed and measured field 
64 0 change are fairly easy to detect and eliminate. 

641 Recall from earlier sections that two length scales were usually required for a{s) — a 

642 shorter one of a few tens of meters or less to fit the peak measured field change, E mp , and a 

643 longer one of a few hundred meters to fit the amplitude of E m {t) at later times. We have tried to 
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64 4 capture the relative importance of these two scales by focusing above on the 100 m and 1000 m 
645 levels, but it can be seen better by examining the relevant fit parameters in Table 1 . Mean values 

64 6 of these parameters over the 1 9 strokes for which all four of them are defined are a mi „i = 0.72 ± 

647 0.09, L a i = 31 ±29 m, a mi „2 = 0.62 ± 0. 15, and L a 2 = 430 ± 280 m. As mentioned previously 

64 8 (e.g., Figure 12), the very rapid decreases in a(s) with increasing path length that are inferred in 

64 9 several cases imply that these return strokes deposit a significant fraction of their total charge in 

650 the lowest few tens of meters above the surface. 

651 Such localized charge depositions are not implied by very close measurements, however. 

652 For example, Schoene et al. [2003] found that the ratio of return-stroke electric-field change 

653 measured 1 5 m from the channel base to that at 3 0 m range averaged 1 .76 ± 0. 1 5 over 77 rocket- 

654 triggered strokes, suggesting a range dependence of r' 0,8 — a fall-off noticeably slower than Hr. 

655 This may be compared with a similar average ratio for the field change due to the immediately 

65 6 preceding dart leaders of 1.88 ± 0.15, suggesting a range dependence of r' 0 9 . Such a slow fall- 

657 off of field change with range could be interpreted to indicate that the magnitude of the linear 

658 charge density deposited by the stroke (usually assumed to be equal and opposite to that 

659 deposited by the leader) increases with height over the lowest hundred meters of so of channel 

6 60 [Rakov et al., 1998]. In contrast, onr own calculations of the close field changes that would be 

6 61 produced by model strokes having very rapid decreases in a(s) (such as 8725/3 and 8726/3) 

662 predict a fall-off with range that is considerably faster than Hr within a few tens of meters of the 

663 channel base. 

664 We have found only one reference in the literature that appears consistent with our 

665 inf erence of rapidly decreasing a(s) near the surface. Wang et al. [1999] reported that time- 

666 resolved optical imaging near the channel base of two triggered strokes showed peak RLI to 

667 attenuate by about 30% over the lowest tens of meters. In spite of this lone observation, 

668 however, the discrepancy between our inference and the recent observations of close electric 

669 fields cited above leads us to question the meaning of our deductions. An alternative 

67 0 interpretation that is motivated by our streak images and appears more physically reasonable is 

671 proposed in the next section. 

672 Further insight into the need for coupling a value of a mi „i that is significantly less than 

673 unity with a very short L a j in certain cases (especially 8715/7, 8715/9, 8715/10, 8725/3, 8725/5, 

674 8726/1, 8726/3, and 8726/5) can be obtained by comparison of the physical propagation speeds 
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675 that are discussed above with the so-called transmission- line-model velocity, vtlm= - 
67 6 (2tiD/ \xo)E mp /i p (0) [e.g . , Willett et al., 1989]. This v T lm is just another way of looking at the ratio 

677 of peak radiation field measured at horizontal range, D, to peak current measured at the surface. 
67 8 Its values have been recalculated at our observing range, D- 5.2 km, using the new zero and 

67 9 peak levels of E m (t) that are shown in Figures 1 5, where we have attempted to better remove any 

68 0 field change due to the leader. The ratio of observed V 2 D to calculated v T lm averages 1 . 17 ± 0.20 

681 over the 1 4 strokes for which both are available. The fact that V 2 D tends to be the larger of the 

682 two indicates that its use as a physical propagation speed for the return-stroke current front in our 

683 (or any other transmission-line-like) model will tend to over-predict E mp in many cases (as 

684 already found, though to a lesser extent, by Thottappillil and Uman [1993]) unless draconian 

685 measures are taken to prevent this. 

68 6 We have already seen in the previous section that the inferred 3-D speed, v e fj(s), when 

687 averaged over most of the reconstructed channel, tends to be even larger than v?/> Thus, it is not 

688 surprising that v e jj( 1 0)/vtlm (where the effective front speed at the 10 m level from Table 2 has 
68 9 been used as the most relevant to the peak of E c (t)) averages 1 .34 ± 0.32 over all 28 return 

6 90 strokes. This implies an even greater over-prediction of E mp in the present modeling, as 
6 91 illustrated in Figure 4 for example. To obviate this over-prediction, a very rapid decrease in a{s) 

6 92 with increasing path length is required in several cases. Indeed, the same strokes listed in the 

693 previous paragraph have among the largest ratios, v e jj( 1 0)/vtlm- 

6 94 Inferred peak currents at several levels have been listed in Table 4 for the 14 strokes also 

6 95 having streak photographs. Mean values of i p (s ) are given in Table 5 for comparison with the 

696 corresponding average amplitudes of relative light intensity from Table 3 (shown in the same 

697 manner as those of 10 - 90% RLI, discussed above). It is gratifying that both of these averages 

698 decrease at roughly the same rate with increasing height. (The agreement becomes nearly 

6 9 9 perfect if the mean RLI amplitude at optical level one is compared with i p at the surface instead 

7 00 of that at 30 m. An argument for this adjustment will also be discussed in the next section.) 

7 01 Looking back at the left sides of Figures 1 8, we also see that individual profiles of i p (s) and RLI 

7 02 amplitude agree reasonably well, both in slope and in overall relative magnitude, although 

7 03 exceptions to the former include strokes 8715/10 and 8726/1 and exceptions to the latter include 

7 04 8726/2, 8726/3, 8726/4, 8732/1, and 8732/2. 
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705 The evolution of relative light intensity with height in return-stroke channels has also 

70 6 been evaluated photographically by Jordan and Uman [1983] and by Jordan et al. [1997]; Mach 

707 and Rust [1989] do not present formal analyses, but assert that their photoelectric data is 

708 generally consistent with that of Jordan and Uman [1983]. Fundamentally, the peak luminous 

7 09 intensity of the subsequent return stroke is observed to decrease with height by about a factor of 

710 two over 500 m, in fair agreement with the present results. 

711 Another noteworthy feature of the left sides of Figures 1 8 is that the RLI amplitude at 

712 optical level one is occasionally much brighter than that at level two. That is, there is sometimes 

713 a dramatic change in the slope of the amplitude profile at level two. This is especially evident in 

714 strokes 8715/9, 8726/2 and 8728/10. Looking at the streak photographs themselves, it is obvious 

715 that the first and last of these three strokes, plus 8725/3 and 8726/1, are preceded by particularly 

716 bright dart leaders. Figure 19 shows the most dramatic example, stroke 8725/3, which may he 

717 compared with more typical stroke 8725/1 in Figure 17. Since the optical emissions from leader 

718 and return stroke are generally indistinguishable at optical level one, because of the limited time 

719 resolution of the streak camera, the tabulated values of stroke RLI amplitude at the lowest 

720 measurement level must have been significantly exaggerated by leader light emission in at least 

721 some cases. 

722 An objective evaluation of any such exaggerated stroke brightness at the lowest optical 

723 level would involve a comparison of the brightness of the leader relative to its return stroke. 

724 This comparison can be made reasonably quantitative by forming the ratio, fiom Table 3, of the 

725 RLI "baseline" to its "amplitude" at each level. Recall that this baseline is measured 

72 6 immediately before the onset of the return stroke and thus approximately represents the residual 

727 brightness of the leader channel (relative to a zero at the same level before leader onset), whereas 

728 the amplitude is the difference between the peak RLI during the onset of the return stroke 
72 9 (relative to the same zero) and the corresponding baseline value. (At the higher levels this 

730 baseline value may be substantially less accurate than the stroke amplitude because considerable 

731 time may have elapsed since the corresponding zero reading. Note that the leader tip is often 

732 appreciably brighter, as in Figure 21, but is not measured here.) This ratio has been calculated 

733 for each stroke at every level (except where "baseline" is negative or where "amplitude" is small 

734 enough to produce excessive noise in the ratio) and then averaged over all levels except the 

735 lowest (where leader and return stroke are presumed indistinguishable). The results are 
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7 36 indeterminate in two cases (8725/5 and 8726/3) because of low return-stroke amplitudes, but 

737 otherwise range from near zero (strokes 8726/2, 8726/4, 8732/1 and 8732/2) to over 30% 

7 38 (8715/9, 8725/3, 8726/1, and 8728/10), the largest ratio being 59% for stroke 8725/3. 

7 39 The first indication that anomalously bright dart leaders do occur was from Guo and 

7 40 Krider [1985], who reported one or more of them in about 5% of natural multiple-stroke flashes 

741 in Florida and argued that leader brightness might occasionally equal or even exceed that of the 

7 42 corresponding subsequent return stroke. Idone and Orville [1985] found that the ratio of dart- 
7 43 leader to return- stroke RLI (measured 50 m above the channel base) averaged 0. 1 ± 0.07 over 22 

7 44 triggered strokes in which leader brightness could be measured; and they estimated the 

7 4 5 corresponding ratio of peak currents to average 0. 17, but with a rather large range of [0.03, 0.3], 

7 4 6 Jordan et al. [1997] reported that the brightest three of 23 dart leaders in natural subsequent 

7 47 strokes had peak RLIs 30 - 50 % that of the corresponding return strokes. They also referred to 

7 4 8 (but did not report) measurements of the "plateau" brightness — after the dart peak but before the 

7 4 9 stroke onset — which is the quantity that is approximated by the "baseline" values in our Table 3. 

7 50 Mach and Rust [1997] reported what they called "postdait" brightness relative to that of the dart 

751 peak. Dividing mean values from their Table 1, we infer an average ratio of postdart to return- 

752 stroke peak brightness of roughly 1/16. Finally, Kodali et al. [2005] reported a mean ratio of 

7 53 dart-leader current (inferred from their near-field measurements on triggered strokes, assuming 
7 54 the leader current above the tip to be uniform and constant) to measured return- stroke peak 
75 5 current of 0.22. Based on these references, the values in our Table 3 seem reasonable. 

756 

7 57 Leader vs. Return-Stroke Current 

758 

759 Table 6 compares our objective measure of leader/retum- stroke relative brightness with 

7 60 two possible indicators of draconian measures to fit E mp — a(100) and v^(10 )/vtlm, both 

7 61 discussed in the previous section. As expected, the correlation coefficient between a(100) and 

7 62 v e ^(10)/vriMis high - -0.88 over all 24 strokes, significantly different from zero at well over the 

7 63 99% confidence level. The correlations of our leader/retum-stroke brightness ratio with these 

7 64 two parameters are -0.74 and +0.41, respectively, over 12 strokes. Thus, the brightness ratio is 
7 65 significantly correlated with the 100 m charge- amplitude multiplier at the 99% confidence level, 
7 66 although its correlation with the velocity ratio is not statistically significant. These results 
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7 67 suggest that the peak currents measured at the surface exaggerate those in the return strokes 
7 68 themselves by essentially the same mechanism described above, at least in cases of relatively 
7 69 bright (high-current) leaders. We suspect that this is the reason for the above-mentioned 
77 0 tendency toward over-prediction of E mp , as well as for the extremely rapid decreases in a(s) that 
771 we have employed to obviate it. 

7 72 Although this is a new idea, to our knowledge, the relationship between leader and 

7 73 return-stroke current that it presumes is not new. For example, Lin et al. [1980] assumed that the 

7 7 4 leader current continues to flow after the onset of the return stroke, remaining uniform and 

7 75 constant above the advancing return-stroke front. This concept is re-iterated by Jordan et al. 

77 6 [1997], who remarked that the optical "plateau" that they observed in natural dart leaders 

111 suggested "that a steady leader current flows through each channel section behind the downward 

77 8 moving leader tip before, and perhaps for some time after, the return stroke has passed that 
7 7 9 channel section." 

7 80 Consider the following simple argument, based on the transmission-line model. (Here we 

781 ignore the transient effects of any upward-connecting discharges and/or reflections at the ground 
7 82 on the currents and radiation fields in these triggered subsequent strokes.) Let the current and 
7 83 front-propagation speed of the dart leader be fr, and -v L , respectively, while the corresponding 

7 8 4 parameters for the return stroke are is and +vs- Near the surface the leader front radiates a field 

7 85 proportional to -frv/, that turns off when the leader reaches ground. The return stroke then begins 

78 6 radiating a field proportional to -isvs as it propagates up the leader channel from ground. 

7 87 Assume, with the references in the previous paragraph, that the leader current continues and 

7 88 remains uniform along the channel after stroke onset. The peak total current at any level 

7 8 9 (including that measured at stroke onset by a shunt at the surface) is then + is, which may be 

7 90 significantly greater than is alone in cases with relatively bright leaders. Since fr < is and vi « 

7 91 v* however, the change in radiated field magnitude at return-stroke onset is proportional to -i s vs 
7 92 - (-frvi) ~ -i s vs- (For example, taking the average ratios, v L ~ vs/6 from Idone et al. [1984] and k 

7 93 ~ 0.22 is from Kodali et al. [2005], we find a radiation-field change proportional to -0.96frv,s'.) 

7 94 Thus, it is essentially only the return-stroke current itself that is relevant to the distant (radiation) 

7 95 field change. Our measured E m (t) at D - 5.2 km consists essentially of radiation field at the time 

7 96 of E mp , so this analysis should be relevant here. 
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7 97 The above argument appears consistent with both the streak photographs and our inferred 

7 98 rapid decreases of a(s ) and i p (s) in cases of relatively bright leaders. We claim that the latter are 

7 99 merely the inevitable result of presenting our model-fitting procedure with exaggerated 

8 00 measurements of retum-stroke-current amplitude at the surface. It even makes sense to compare 

801 the mean RLI amplitude at optical level one with i p at the surface if both measurements then 

802 fully include the effects of the leaders. If correct, this argument obviously also has implications 

803 for the widespread estimation of return-stroke speeds from measured E mp and i p ( 0) through v T lm- 

804 This matter could be further investigated using modem, near-field measurements, such as have 

8 05 been reported recently from Camp Blanding, Florida [e.g., Kodali et al., 2005], together with far- 

8 06 field measurements of the peak fields radiated by the same return strokes. 

807 The question remains whether the above inferences about a(s) and i p (s) can be corrected 

808 for this effect. Since we do not know how to accurately remove the leader current from the 

8 0 9 measured surface-current waveforms, however, we confine ourselves here to a single illustrative 

810 example. In stroke 8725/3 both the velocity ratio, v e jj(l0)/vTLM = 1 -62, and the level- averaged 

811 leader/retum-stroke brightness ratio, 0.59, suggest that the peak current due to the return stroke 

812 alone should be reduced to about 62% of the measured peak current at the surface. In the model, 

813 a slightly larger reduction of this peak current to 56% of the measured value allows us to entirely 

814 elim in ate the shorter of the two length scales for a(s). (To achieve this reduction, the measured 

815 waveform is decreased to 56% of its original magnitude throughout the time interval before its 

816 peak, whereas a constant 18.9 kA is subtracted from it at all later times. Obviously this method 

817 of modifying the current waveform is somewhat arbitrary.) A fit of E c (t) to E m (t) that is 

818 essentially the same as that shown in Figure 15c can then be obtained using a min i = 0.82, L a j = 

819 500 m, and a mi „2 =1.0 (all other parameters remaining as in Table 1). The resulting 

820 measurements of 10 - 90% current rise time and of effective stroke-front speed at the various 

821 levels remain essentially unchanged from those in Table 2. Although a(s) changes dramatically, 

822 of course, i p (s) changes from the values in Table 3 only at the lowest two levels, remaining 

823 essentially the same at 30, 100, 300, and 1000 m. This result is considered to corroborate our 

824 suspicion of exaggerated return-stroke currents near the surface. 

825 Thus, the inference in Tables 1, 2, 4, and 5 (and elsewhere) of very rapid decreases in 

82 6 ( 3 ( 5 ) and i p (s), and the consequent concentrations of charge deposition that are illustrated in 

827 Figures 7 and 12, near the surface in several cases are probably invalid. These artifacts are 
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828 argued to result from using the total measured current at the surface as the initial current into the 

82 9 base of the return stroke itself. There is no fully satisfactory way to correct this problem in the 

830 present dataset. Nevertheless, it appears that the inferred values of current rise time at all levels 

831 and of average effective stroke-propagation speed are unaffected, even in strokes with very 

832 bright leaders. Further, we believe our inference of decreasing current amplitude above about 

833 100 m of path length (taken as a compromise between the mean values of L a i -31m and L c ,2 = 

834 430 m, noted in the previous sub-section, the former length scale now being assumed to be 

835 largely spurious). The last two rows of Table 4 are intended to better distinguish the presumed 

836 valid and the probably invalid variations of inferred peak current. There we have first scaled the 

837 values tabulated for each stroke relative to that at 1 00 m and then calculated means and standard 

838 deviations of these scaled results for each level. We conclude that the peak current probably 

839 does decrease by about 37%, on averaged, between 100 m and 1000 m of path length above the 

8 4 0 ground. On the other hand, the peak current may well increase much less than 58%, on average, 

841 between 1 00 m and the surface. 

842 

843 

8 44 Summary and Conclusions 
845 

84 6 This paper has attempted to take advantage of three-dimensional channel reconstructions 

84 7 of rocket-triggered lightning flashes in Florida in order to infer the behavior of the current in 
8 4 8 return strokes above the ground from current waveforms measured at the channel base and 

84 9 electric-field-change waveforms measured at a range of 5.2 km. The reconstruction of six 

8 50 lightning channels from stereo photographs was described in detail in Appendix A. The method 

851 of calculating field-change signatures due to return strokes that follow these piecewise-linear 

852 channels was then explained. Formal assumptions about the variation stroke-current waveforms 

853 along the channels constrained the problem enough that the rise time, propagation speed, and 

8 54 peak amplitude of the current could be estimated as a function of path length with reasonable 

855 confidence. Results were first presented on the tacit assumption that the measured current 

8 5 6 waveforms were entirely due to the return strokes themselves. 

857 Some of the results derived in this way that we still believe are the following: 1) The fine 

858 structure of the field-change waveforms that are radiated by subsequent return strokes can be 
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85 9 explained, in large part, by channel geometry, although it can also be affected significantly by 

8 60 the shape of the current waveform that enters the channel base. 2) The average 10 - 90% rise 

8 61 time of the stroke current increased by about a factor of seven in our sample of 24 triggered 

8 62 strokes, from an observed 0.3 1 ± 0.17 ps at the surface to an inferred 2.2 ± 0.5 ps at 1 km path 

8 63 length above the surface. 3) The three-dimensional propagation speed of the current front 

o 

8 64 averaged 1 .80 ± 0.24 X 1 0 m/s over channel lengths typically greater than 1 km for the same 24 
8 65 strokes. 

8 66 Next, streak photographs of a subset of these return strokes were analyzed in terms of 

8 67 relative light intensity versus path length and time, and the resulting estimates of rise time, 

8 68 propagation speed, and peak amplitude of the. RLI were compared with the corresponding 
8 69 electrical estimates for current. Although these comparisons appeared generally reasonable, two 
8 70 anomalies were noted that suggested a variation in the modeling. First, it was remarked that 
8 71 even the optically measured, two-dimensional, stroke-propagation speeds (which were an 
8 72 appropriate factor smaller than the electrically inferred, three-dimensional, current- front- 
873 propagation speeds) led to an over-prediction of peak radiation fields in several cases, unless 

87 4 draconian measures (which appeared to conflict with other results in the literature) were 

875 employed to obviate this. Second, it was noted that peak RLI amplitudes for return strokes 
87 6 sometimes decreased dramatically between the lowest measurement level, where the leader and 
877 return stroke could not be distinguished, and the next higher level, where they could be. Both of 
8 7 8 these anomalies tended to occur most markedly in cases where the leaders were brightest relative 

8 7 9 to their strokes. It was concluded that (4) a significant fraction of the measured current at the 

8 80 surface was probably due to the leaders in cases when they were relatively bright. Therefore, our 

8 81 assumption that these measured currents were entirely due to the return strokes was forcing an 

882 unreasonably large and abrupt reduction in inferred current amplitude over the first few tens of 

883 meters above the surface. 

8 84 With this conclusion in mind, the first inferences of current amplitude as a function of 

8 85 height were re-examined. It was judged that the anomalously abrupt decreases with increasing 

886 height near the surface were probably spurious, but that the slower increases at higher altitudes 

887 were probably still valid. Thus, the final conclusion of our study was that (5) return-stroke peak 

8 8 8 currents decreased by about 37 ± 12% from 100 m to 1 km of path length above the surface over 
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88 9 our 24 strokes. Although these peak currents likely also decreased between the surface and 1 00 

8 90 m, this decrease was probably not as great as originally inferred. 

891 

892 

8 93 Appendix A — 3-D Channel Reconstruction 
894 

8 95 Figure A1 shows one example (flash 8732) of the six stereo pairs of still photographs that 

8 96 were taken with automated 35 mm cameras from the two sites (NRL and SUNY) during the 1987 
8 97 experiment. These photographs were all printed at the same magnification (6.8X) on 8" X 10" 

8 98 paper, and the lightning channels were then digitized by marking a point at every visually 

899 detectable deviation from straight lines. Typically, this produced a few hundred points per 

900 channel image that were tabulated in centimeters ( X , Y) relative to the "principal point" [Hallert, 

901 1960, p.18], which was taken to be at the center of the print (±0.5 cm). Connecting these points 

902 with straight lines results in a piecewise-linear representation of the channel. 

903 Figure A2a shows the locations of the two cameras relative to the triggering site, which is 

904 taken as the origin of a right-handed, Cartesian coordinate system (x, y, z ) with its z-axis vertical 

905 and its x-axis passing through the NRL site (almost due east of the triggering site). The position 

906 parameters for each site are given in the figure, as determined (±0.6 m) from a survey of the area 

907 that was conducted by KSC. The azimuthal orientation of the cameras (rotation about the 

908 vertical axis passing through each camera — e.g., ys in Figure A2b) could not be accurately 

90 9 measured in the field, but the "swing" (rotation of the image frame about the "optical axis" — the 

910 ray passing through the optical center of the lens and the principal point on the image) was 

911 approximately zeroed (±0.3 deg.) by leveling the camera bodies, and the camera elevation angles 

912 (rotation about the horizontal axis that is perpendicular to the optical axis — e.g. , o. N in Figure 

913 A2c) were estimated (±0.5 deg.) with a spirit level to be 10.0 deg. at NRL and 9.5 deg. at SUNY. 

914 Finally, the focal lengths of the camera lenses (nominally 55 mm at NRL and either 35 mm or 24 

915 mm at SUNY) were carefully measured (±1%) to be 54.7 mm, 35.7 mm, and 23.7 mm, 

916 respectively. The "barrel distortion" of the 24 mm lens that produced all of the SUNY 

917 photographs except that of flash 8715 was also measured and was applied to correct the 

918 corresponding tabulations (but see the discussion below). 
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919 Extrapolating the lightning channels (which were nearly straight at the bottom, following 

92 0 the triggering wires) to ground level (estimated from the measured camera elevation angles), it 

921 was possible to accurately determine the azimuth angles of the two camera optical axes ( e.g . , ys 

922 in Figure A2b). This procedure was insensitive to small errors in the camera elevation angles 

923 because the triggering wires were nearly vertical. By visually identifying the top of the 

924 triggering wire on each image and forcing these points to coincide in the 3-D reconstruction of 

925 each channel, it was also possible to correct the measured elevation angle of the NRL camera (a,y 

92 6 in Figure A2c — note that this sketch has been slightly simplified and gives the precise definition 

927 of cxn only when y,v = 0 so that the optical axis lies in the x-z plane) relative to that of the SUNY 

928 camera, which was assumed correct. It would have been better to have two or more fixed objects 

92 9 of kn own location in the field of view of both cameras, from which to fully determine the camera 

930 orientations, but the triggering site was not suitable for this refinement. Table A1 lists the 

931 camera-orientation parameters that were determined in this manner. Notice that the deduced 

932 camera azimuths, as well as the corrected elevation of the NRL camera, varied slightly from 

933 flash to flash, presumably due to instability of the tripods on which they were mounted. 

934 The formulae of Hallert [ 1 960, Appendix A], together with the camera-orientation 

935 parameters tabulated above, have been used to convert the points (X, Y) that were digitized from 

936 the photographic prints into absolute angle pairs (0, cp) — elevation measured upward from 

937 horizontal and azimuth measured clockwise from the direction of the y axis, respectively — that 

938 define rays from the camera positions. Then the position parameters in Figure A2a have been 

939 used to compute the intersections of these rays (x, y, z) in our Cartesian coordinate system. The 

940 latter step is non-trivial and was performed in practice by solving simultaneous equations with 

941 the Mathematica (TM) software, following an outline originally due to Stan Heckman [personal 

942 communication, 1995]. 

943 Figure A3 is a conceptual illustration of the channel-reconstruction procedure. (Note that 

944 a P here has been given the opposite sign from Figure A2a, solely for convenience of illustration.) 

945 The central problem is that, in general, it is not possible to identify any given point on the 

94 6 ligh tnin g channel in one image with the corresponding point in the other image. Therefore, this 

947 is not the classical double-theodolite problem, where there is redundant data (four angles to 

948 determine a point in 3-space). A heuristic description of the channel-reconstruction method is as 

94 9 follows. The ray, /,, that is defined by point r in the first camera image can be projected onto the 
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950 second camera image as the line, which intersects the second channel image somewhere (in 

951 general, between two of the vertices, ,v and j+1, as shown in the inset of Figure A3). This 

952 intersection point — call it s' — defines a ray from the second camera, l s >, that intersects the first 

953 ray, defining a point in 3-space — identified by the indices, (r, s') — on the reconstructed channel. 

954 The apparently redundant piece of information (the fourth angle) is effectively used up to 

955 determine the intersection point, s', on the second channel image. 

956 Mathematically, our solution of this problem proceeds from two vector equations, each 

957 specifying the location in 3-space of the intersection point, r, of rays l r and l s < in terms of the 

958 angles, (0,, (pi), and the range, n, from camera i: 

959 

960 

r = r<)i + ri[xsin (0i) cos (0i) + y cos (0i) cos (0i) + z sin (©i ) J 

(Al) 

961 r = r 02 + r 2 [x sin (0 2 ) cos (0 2 ) + y cos (0 2 ) cos (0 2 ) + z sin (0 2 ) ] 

962 

963 

964 The #*&• are the positions of the two cameras; x , y , and z are the unit vectors in our Cartesian 

9 65 coordinate system. Equating these two expressions for r, we obtain three linear scalar equations 

9 66 in two explicit unknowns, r } and r 2 . Recall, however, that the angles from camera 1, (0j, (pi), are 

9 67 implicitly functions of the index, r, whereas those from camera 2 are implicitly functions of s — 

968 we have linearly interpolated the angles between integral values of these indices — and that it has 

969 not yet been explained how to find the s' (see Figure A3) that corresponds to any particular value 

970 of r. Given r, we might simply solve numerically the three simultaneous equations (now non- 

971 linear in the implicit variable, 5 ) for n, r 2 , and s', as Stan Heckman [personal communication, 

972 1995] originally suggested. Once rj and r 2 are known, it is obviously trivial to find r from either 

973 of Equations Al . (This approach might be generalized to find the "best" solution, in the least- 

97 4 squares sense, from three or more camera images of a single lightning channel.) Because we 

975 have only two cameras, however, the computation can be simplified as follows: 

97 6 Three linear equations in two unknowns have a unique solution if and only if the rank of 

97 7 the 3X2 coefficient matrix and the rank of the 3X3 "augmented matrix," A (containing an 
978 additional column of the constant terms), both equal two [e.g. , Boas, 1 966, Section 3.7]. Thus, 
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the determinant of the augmented matrix must equal zero. Given a particular value of r, the non- 
linear equation 

det(y4) = 0 (A2) 

can be solved numerically for the implicit variable, s'(r). From Equations A1 the augmented 
matrix is found to be 


A = 


/-sin ( 0 1 ) cos (0i) sin (0 2 ) cos (0 2 ) x 0 i - x 02 

-COS ( 01 ) COS (Si) COS ( 0 2 ) COS (0 2 ) Yoi - y02 

-sin (0i) sin (0 2 ) z 0 i - z 0 2 ) 


(A3) 


where (xoi,yoi, z od are the known Cartesian components of rou the angles, (6f, <p } ), are known 
functions of the specified index r, and (0 2 , <p 2 ) are known functions of the unknown index, s. 
Once A is determined, the two linear equations with the largest determinant are solved for r 2 and 
r 2 , and r is foxmd. 

Occasionally, Equation A2 has multiple solutions for s' at a given r. One way to deal 
with this situation is as follows. Notice that the two rays, one from each camera intersecting at 
the reconstructed channel point (r, s') in Figure A3, define a plane. This "solution plane" is fully 
determined by the locations of the two cameras and by the point, r, on the first channel image, as 
described heuristically above. The plane cuts each camera image in a projection line, as further 
illustrated by l' r and l' s > in Figures A4, where these projection lines correspond to the point, r = 
155, on the NRL channel image of flash 8732. [Channel points are numbered here in ascending 
order from the channel base upward. Note that Figures A4 actually show the "image planes" 
after transformation from the original linear measurements on each print, (X, Y), into the absolute 
angles (0, <p) — independent of camera orientation, focal length, etc. — since this was the notation 
in which the data were analyzed mathematically to find points on the lightning channel. Thus, 
"lines" l'r and r s < in Figures A4 are not actually straight, although they appear so in these small 
sections of the images.] There are two interesting facts to note here: (1) Any channel-image 
vertices that lie on the projection line in one photograph project onto the same line in the other 
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1007 photograph and are thus indistinguishable. For example, points 152 - 156 on the NRL image all 

1008 correspond to the projection line, shown on the SUNY image of the same flash, which 

100 9 intersects that channel image near the point, s = 145, among other places. Close examination of 

1010 the original photographs suggests channel propagation almost directly toward the SUNY camera 

1011 in this region. (2) Downward apparent propagation relative to the projection line in one camera 

1012 image must correspond to downward apparent propagation relative to the corresponding 

1013 projection line in the other image. Hence, certain channel kinks can be unambiguously identified 

1014 in both images. For example, the downward loop that is defined by points 162 - 176 in the NRL 

1015 image corresponds to the downward loop, points 150 - 163, in the SUNY image. 

1016 Using these conclusions, we can resolve the apparent ambiguity of the multiple solutions 

1017 that are indicated by the red [projection lines in Figures A4. Points 152 - 156 on the NRL image 

1018 correspond to points 145 and 146 in the SUNY image (not to point 149 nor point 164), whereas 

1019 NRL segment 161 - 162 corresponds to SUNY point 149, and NRL segment 176-177 

1020 corresponds to SUNY point 1 64. 

1021 Loops of the same sense, such as the pair (NRL 162 - 176 and SUNY 150 - 163) 

1022 identified above, that could be unambiguously associated between the two images were used as a 

1023 check on the uncertainties that are inherent in our 3-D channel reconstructions. In most of the 

1024 flashes the top of the triggering wire was readily identifiable in both photographs. (As 

1025 mentioned above, the elevation angle of the NRL camera was adjusted — in each case by less 

1026 than the 0.5 degree uncertainty in our measurement of that angle — to make the top of the wire in 

1027 both images coincide in the reconstruction.) In each of two flashes — 87 17 and 8732 — there 

102 8 were two obvious ki nk s, in addition to the top of the wire, that could be used to check the 

102 9 reconstructions. (Figures A4 illustrate one of these four major kinks.) Surprisingly, best 

1030 agreement was obtained by eliminating the correction for barrel distortion of the 24 mm lens at 

1031 the SUNY site. In fact, it was found that all identifiable kinks in all of the channel images could 

1032 be made to coincide by this simple parameter change. In retrospect this seems reasonable, since 

1033 the barrel distortion was measured at close focus, whereas the lens was actually used at infinity 

1034 focus, where aberrations are usually minimized by design. Therefore, it is presumed that the 

1035 reconstructions without the barrel-distortion correction are the best possible under the 

103 6 circumstances. (The data in Table A1 correspond to this assumption.) 
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1037 The significance of eliminating the barrel-distortion correction, and an example of the 

1038 uncertainty that is inherent in our reconstructions, are illustrated by Figures A5a & A5b, which 

1039 show two views of reconstructed channel 8732, with (blue) and without (red — preferred) the 

1040 correction. Notice that the differences between these two reconstructions are a few tens of 

1041 meters or less throughout, which is typical. Two views of all six reconstructed channels are 

1042 given in Figures A6a & A6b. Flashes 8715, 8717, 8725, 8726, 8728, and 8732 are shown in 

1043 cyan, magenta, blue, green, red, and black, respectively. Notice that there is considerable 

1044 variability - much greater than a few tens of meters - in both the length of the triggering wire 

1045 and the overall channel shape among these flashes. 

10 4 6 In order to illustrate the impact of uncertainties in the reconstructed channels, Figure A7 

1047 compares electric-field changes for stroke 2 of flash 8732 that have been calculated from the two 

1048 different channel reconstructions shown in Figures A5. The shape, amplitude, and propagation 

10 4 9 speed of the return-stroke current waveform that is deduced for this event in the body of the 

1050 paper has been used in both cases. Again, the blue curve is with, and the red (preferred) curve is 

1051 without, the barrel-distortion correction. The differences between these two waveforms are 

1052 small and are typical of the impact of geometrical uncertainties on the model field changes in our 

1053 dataset. This satisfying result is a consequence of both the good overall accuracy of the 3-D 

1054 reconstructions and the considerable smoothing of fine structure that is caused by the rapid 

1055 increase in current rise time with height, as discussed in the body of the paper. 

1056 
10 5-7 

1058 Appendix B — Calculation Method for Piecewise-Linear Channel 

1059 

10 60 A numerical code was written in Mathematica [TM], following the general outline of the 

10 61 FORTRAN code that had been developed previously by Le Vine and Meneghini [1978], to 
10 62 compute the electric-field change at the observing site from simulated return strokes in the 
10 63 reconstructed channels. The adoption of the "generalized TLM" form for the stroke current 
10 64 (basically, Equations 5-7) allowed a number of simplifications in these calculations. The key 

10 65 attributes of this model are that the entire current waveform propagates monotonically upward 

10 66 along the channel while its principal parameters, a(s), v(s), and K(t', s ), depend only on position. 

10 67 Thus, the "pure-TLM" current parameters that are required on each linear channel segment for 
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1068 Equations 2-4 (the fixed current wave shape, IrLtAfirgiiment}, and constant propagation speed, 
10 69 v) can be computed in advance from knowledge of the channel geometry. Then the field-change 
107 0 calculation can proceed time step by time step. 

10 71 At each successive time step the current integral on each channel segment is updated for 

1072 use in Equation 4, and a list of the segments from which radiation can reach the observer is 

1073 computed. Then the contributions to the total field change are summed, both over Equations 2 - 

1074 4 arid over all such "radiating" segments. In this way the field-change waveform is built up over 

1075 time. The size of the time steps is not critical (as long as they are short enough to compute the 
107 6 current integrals with sufficient accuracy), except in the sense that waveform details will be 

1077 missed if they occur entirely between time steps. As long as the channel segments are made 

1078 short enough, this calculation method has been shown to accurately approximate the exact field 

107 9 change from a tortuous channel, even when the peak amplitude, rise time, and propagation speed 

1080 of the current waveform all depend strongly on position. The code has even been shown 

1081 accurate in the near field by comparison with independent numerical calculations of Jens 

1082 Schoene [personal communication, 2003] of the University of Florida. 

1083 

1084 

1085 Appendix C — Analysis of Streak Photographs 

1086 

10 87 Relative light intensity (RLI) determinations from the streak recordings were possible 

10 88 because the data strips were developed with a "calibration strip" of known relative exposure. 

10 8 9 Each calibration strip was exposed to a single Xenon flash of about 3 us duration, which evenly 

10 90 illuminated a film strip positioned directly behind a Kodak calibration step tablet. The tablet has 

10 91 21 steps of known density values, thereby transmitting a lOOOx range of RLI. With this 

1092 information, film density can be converted reliably to RLI on a microsecond time scale, as 
10 93 previously demonstrated convincingly by Jordan and Uman [1983], 

10 94 Here, a Xillix 1412 CCD camera was used to image the data and calibration film strips. 

10 95 The Xillix camera has an image-plane sensor of 1344 x 1035 square pixels and 12-bit output, 

10 96 with a specified dynamic range of >60 dB. Various tests confirmed that a factor of slightly more 

10 97 than a thousand in RLI could be reliably recognized with this device, comparable to the range of 

10 98 the Kodak step tablet and the film emulsion itself. The data and calibration strips were imaged 
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10 99 with the Xillix camera under identical illumination conditions; the Xillix output for known 

1100 relative exposures on the calibration strip establishes a "lookup table" for values on the data strip, 

1101 thereby yielding RLI values for the streak image data. Fortunately, the densities for almost all 

1102 the data and calibration strips were outside the problematic "toe" and "shoulder" portions of the 

1103 film-response curve, allowing good interpolation accuracy. 

1104 All images rendered with the Xillix camera were scanned at 77 pixels per mm on the 

1105 film, appropriate to having the film grain ultimately limit the analysis. This yields a vertical 

1106 (spatial) scale of about 0.84 pixels per meter in the object plane and a typical temporal scale of 

1107 about seven pixels per microsecond for the streak camera. 

1108 In this analysis, eight or nine separate vertical levels were first selected for each image, 

1109 being careful to avoid intense "streak lines" or scratches in the emulsion. Figure 17 for stroke 

1110 8725/1 is a typical example. At each pre-selected vertical level, a horizontal strip of 1344x7 

1111 pixels centered on that level was then extracted from the digital image and averaged vertically 

1112 across the seven adjacent pixels (typically about 8 m of channel height) to reduce grain noise. 

1113 Using the streak camera’s known writing rate and the calibration information, the averaged 

1114 values at each level were converted to a time series of RLI. 

1115 The background RLI before the onset of any perceptible leader illumination was 

1116 estimated separately for the time series of each stroke at each level and used as the zero of RLI 

1117 for that particular time series. A typical example of the resulting "raw" time series is given in 

1118 Figure Cl, corresponding to level seven (counting upward from the lowest level, always 

1119 considered level one) of stroke 8725/1 . This procedure should allow the RLI to be compared 

1120 between different levels of the same stroke and also between different strokes in the same flash. 

1121 (Comparison between strokes in different flashes is potentially somewhat problematic, since they 

1122 are normally on different pieces of film developed on different days, although in principle the 

1123 film calibration should permit such a comparison as well.) 

1124 To determine the true altitude and the path length along the channel that corresponds to 

1125 any given level on a streak image, we projected the corresponding three-dimensional channel 

112 6 reconstruction from the perspective of the streak camera and then visually lined up features in 

1127 the two images. (In many cases it was easier to do this with the leader image when it was 

112 8 visible, since it was often sharper than that of the return stroke.) In this way, for example, level 

112 9 seven of stroke 8725/1 (see Figure 17) was determined to be about 728 m above the surface, or 
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1130 about 886 m of path length of from the channel base. Interpolation was required on the nearly 

1131 strait sections of channel where the lightning had followed the triggering wire, of course. (The 

1132 actual bottom of the channel was obscured by an intervening tree line. Consequently, level one 

1133 was taken as 30 ± 15 m altitude in each case, based on several lines of evidence.) This procedure 

1134 allowed the path length corresponding to each level to be determined to an accuracy of ±50 m or 

1135 better. Based on the average 3-D stroke-propagation speeds estimated in Table 2, this 

113 6 corresponds to a temporal uncertainty of ±0.3 ps or better for these events — certainly adequate 

1137 for our purposes. 

1138 Determination of the 10 - 90% optical rise time and the peak optical amplitude are related 

1139 by the choice of onset and peak RLI for the return stroke under consideration. Both 

1140 determinations were made difficult by the noise on the raw RLI data (e.g.. Figure Cl, already 

1141 vertically averaged as indicated above), which became a greater problem as the stroke amplitude 

1142 decreased. In an effort to minimize this noise without unduly broadening the observed rise time, 

1143 the data were temporally smoothed by the application of a weighted moving average. The 

1144 weighting function was a "cosine bell" [e.g., Willett et al., 1990, Equation A8] with an adjustable 

1145 full width at half maximum (FWHM), normalized to have unity area. It was found 

114 6 experimentally that, as this FWHM was increased for a given time series, the ratio of apparent 

1147 rise time to FWHM approached 0.80 (when the rise time became dominated by the width of the 

1148 weighting function). Therefore, FWHM was kept as small as possible, consistent with the 

114 9 unambiguous determination of rise time, and an absolute minimum of 1 .6 was imposed on this 

1150 ratio. 

1151 In choosing the onset and peak RLI for each stroke, we focused only on the fast-rising 

1152 portion of the time series (e.g., about 75.5 to 82.5 ps in Figure Cl). This was done both to avoid 

1153 broadening of the rise time by any light scattering either in the camera or in the film itself and to 

1154 ignore the gradual rise or "hump" that is often present later in return-stroke light emissions, 

1155 particularly above the surface (again see Figure Cl). In general, these levels were chosen by 

115 6 applying excessive smoothing to find an average "baseline" level just before the onset of the fast- 

1157 rising portion, and an average peak or plateau just after it, that were reasonably independent of 

1158 noise spikes. In the case of large, strongly peaked RLI records (typically at measurement levels 
115 9 near the surface for strokes with large peak currents), however, the peak value was allowed to 
1160 "float" to the maximum of the smoothed waveform as FWHM was varied. 
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1161 Once the onset and peak RLI were determined, 1 0% and 90% levels were "drawn" on the 

1162 smoothed waveform for automatic determination of the last time that the 10% level was 

1163 exceeded and the first time that the 90% level was exceeded, their difference being the estimated 

1164 10 - 90% rise time. The FWHM of the smoothing was then gradually increased from zero until 

1165 the observed rise time stabilized, but not so much that it became steadily increasing. As FWHM 

1166 increased, rise time typically passed through a minimum, which was generally taken as the best 

1167 value, and then slowly increased. Often, however, inconveniently located noise spikes were 

1168 large enough that they caused abrupt jumps in rise time as increasing smoothing caused them to 

1169 fall below one or the other RLI threshold. In these cases the determination of rise time became 

117 0 more subjective. Sometimes two or more values of rise time seemed consistent with the data, as 

1171 illustrated in Figures C2a and C2b, for which the relevant parameters are given in the caption. 

1172 Fortunately, there was normally a rather small range of deduced rise times in such cases. 

117 3 The optical peak amplitude of a stroke was taken as the difference between the baseline 

1174 and peak values determined above. This measurement is also somewhat subjective, the more so 

1175 as the stroke amplitude decreases toward the noise level. 

1176 
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1417 Figure Captions 

1418 

1419 Figure 1 — An example of our "convolution smoothing" of the measured current waveform 

1420 (black curve) to produce different wave shapes, I(ti,s) in Equation 7, at different heights. The 

1421 corresponding values of kernel time scale, r(.v), and 1 0 - 90% rise time are tabulated in the figure. 

1422 

1423 Figure 2 — Two examples of truncated, smoothed, and extrapolated channel-base-current 

1424 waveforms, (a) Full 15 ps time series, illustrating extrapolations, (b) First 1 ps, illustrating 

1425 onset truncation and noise smoothing. 

1426 

1427 Figure 3 — Measured field change for stroke 8732/2 (black) compared to that calculated when 

1428 the measured current waveform propagates at a constant speed of 1.71 X 10 8 m/s without 
142 9 changing either shape or amplitude (green). 

1430 

1 4 31 Figure 4 — Similar to Figure 3, but this time convolution smoothing (r max = 4.36 ps, L x = 7084 m) 

14 32 is included to produce the green curve and then decreasing TLM velocity (v min = 2.1 X 10 m/s, 

14 33 v max = 2.6 X 1 0 8 , L v] = 500 m) is added to produce the red curve. 

1434 

14 35 Figure 5 - TLM-velocity profiles [v(s), solid curves] and current-front effective-velocity profiles 

1436 [v e fj{s), dots] for the two calculations shown in Figure 4 (colors correspond). 

1437 

14 38 Figure 6 — Similar to Figure 4 except that the red curve has been further modified by a 

1439 decreasing amplitude factor (a mi „i = 0.7, L a i = 60 m, a mi „2 = 0.6, L a 2 = 400 m). 

1440 

14 41 Figure 7 - Profile of the amplitude factor [a(s), red curve] and the corresponding relative profile 

14 42 of charge deposition on the channel by the model return stroke shown in Figure 6 (green). 

1443 
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1443 Figure 8 — TLM- and front- velocity profiles, as in Figure 5, for stroke 8732/1 , which required 
14 44 propagation-speed adjustment, as well as amplitude adjustment, in order to match both the peak 
1445 and the subsequent structure of E rn (t). The solid green curve is for v min = 2.0X1 0 8 m/s, v max = 

144 6 2.3 X 10 s , and L v i = 500 m. The red curve adds the paremeters, v amp i = 0.7 X 10 s m/s and L V 2 = 

1447 60 m. 

1448 

14 4 9 Figure 9 — Computed field-change waveforms (colors) corresponding to the velocity profiles in 
14 50 Figure 8 are compared with observed field change (black). 

1451 

14 52 Figure 10 — TLM-velocity (solid) and front- velocity (dots) profiles for stroke 8728/11. The 

14 53 parameters for the solid curve are v min = 2.7 X 10 m/s, v max = 3.0 X 10 , and L vI = 500 m. 

1454 

1455 Figure 1 1 — The effect of a very rapid amplitude-factor decrease in a case with the parameters, 

145 6 t max = 1.16 (IS, L x = 1851 m, v mm = 2.3 X 10 8 m/s, v max = 3.0 X 1 0 8 , L vl = 500 m, and as given in 

14 57 the caption of Figure 12. 

1458 

14 5 9 Figure 12 — Similar to Figure 7 but shown on an expanded height scale for the parameters, a mm! 

14 60 = 0.69, Lai = 5 m, a min2 = 0.35, L a2 = 1000 m, that correspond to the red curve in Figure 1 1 . 

1461 

14 62 Figure 13 — Two different current-waveform extrapolations for stroke 8725/4. 

1463 

14 64 Figure 14 — Model field-change waveforms (green and red, corresponding to the similarly 
14 65 colored current waveforms in Figure 13) are compared with the observed field change (black). 

1466 

14 67 Figure 15 - Observed (black) and modeled (red) field-change waveforms for all 24 return 
14 68 strokes in our data set. (a) collects the events from three flashes for which we could analyze only 

14 69 two strokes each, (b), (c), and (d) each show several strokes from a single flash. The scales are 

1470 V/m versus ps. 

1471 

1472 Figure 16 — Measured and extrapolated current waveforms corresponding to the field-change 
14 73 waveforms of Figure 15. The scales are kA versus ps. 
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1474 


1475 Figure 17 — A typical streak image of leader/retum-stroke sequence for stroke 8725/1 is shown 

147 6 as rendered by the Xillix camera. The analysis levels chosen for this event are illustrated by 

1477 white horizontal bars. The vertical separation between level one (the lowest) and level nine (the 

1478 highest) in this image corresponds to about 980 m at the range of the lightning channel. The 

147 9 horizontal extent of the image corresponds to about 125 ps of time. See Appendix C for further 

14 80 details. 

1481 

14 82 Figure 18 — Comparison of optical and electrical results as a function of path length. The left- 
14 8 3 hand panels show peak-relative-light-intensity and peak-current profiles, while the right-hand 

14 84 panels show 10 - 90% rise-time profiles of RLI and of current. Each flash has its own row of 
14 85 two panels showing all strokes for which both optical data and electrical inferences are available 
14 86 in that flash. Optical results are plotted with solid symbols and lines, whereas electrical results 
14 87 have hollow symbols connected by dashed lines. 

1488 

14 8 9 Figure 19 — Similar to Figure 17 but for stroke 8725/3, which has an especially bright leader. 
1490 

14 91 Figure A1 — Still-camera images of flash 8732 from the NRL site (a), located about 5.2 km east 

14 92 of the triggering site and using a 55 mm lens, and from the SUNY site (b), located about 2.2 km 

14 93 SSE of the triggering site and using a 24 mm lens. 

1494 

14 95 Figure A2 — Camera locations and orientations at the NRL (subscript N) and SUNY (subscript 

14 96 S) sites relative to the triggering site — the origin, O, of the (x, y, z) coordinate system, (a) A and 

14 97 hi represent the relative ranges and heights of the two cameras. The NRL site is taken to be just 

14 98 above the x-axis. Angle ap represents the co-azimuth of the SUNY site, (b) The definition of 

14 99 the camera-azimuth angle at the SUNY site, ys- (c) The (simplified — see text) definition of the 

1500 camera-elevation angle at the NRL site, a n- 

1501 
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1501 Figure A3 — A heuristic illustration of the channel-reconstruction method used in this paper. 

1502 The locations of the two cameras are indicated by Pi (representing NRL) and P 2 (representing 

1503 SUNY, except for a change in the sign of op for convenience of illustration), with the other 

1504 parameters as in Figure A2a. Point r on the channel in the first idealized image defines the ray, 

1505 l r , which projects onto the second image as the line, l' r . This line intersects that imaged channel 

1506 at point, s', which in turn defines a second ray. Ip, and a point on the 3-D channel, identified as 

1507 (r, s') — see text. 

1508 

150 9 Figure A4 — Magnified portions of idealized images from the NRL (a) and SUNY (b) cameras, 

1510 showing a major channel kink, leading to multiple solutions. The graph axes give the absolute 

1511 elevation and azimuth angles, (0„ (pi), in radians. Digitized points on the two images are 

1512 . numbered consecutively from the channel base upward. The "lines", l’ r and l’ s , where the 

1513 "solution plane" corresponding to point r = 155 on the NRL image cuts both images are indicated 

1514 in red. See text for interpretation. 

1515 

1516 Figure A5 - Two views of the 3-D reconstructed channel for flash 8732, with (blue) and without 

1517 (red) the barrel-distortion correction to the SUNY lens. The scales of the (x, y, z) axes are in 

1518 meters relative to the location of the triggering site. 

1519 

1520 Figure A6 - The same two views as in Figure 5, but showing all six reconstructed flashes used 

1521 in this study. 

1522 

1523 Figure A7 - Electric-field changes computed for stroke 8732/2 using the current model inferred 

1524 in this paper (see the section, "Detailed Example of Fitting Procedure") together with the two 

1525 channel reconstructions shown in Figure A5. The blue curve uses the barrel-distortion 
152 6 correction, whereas the red curve does not. 

1527. 

1528 Figure Cl - "Raw" relative light intensity, after averaging seven pixels vertically and after zero 

152 9 subtraction but before any temporal smoothing, at level seven of stroke 8725/1 (see Figure 17). 

1530 
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1531 Figure C2 — Similar to Figure Cl, except that temporal smoothing has been added, the 

1532 "baseline" and peak magnitudes of the fast-rising portion of the return-stroke RLI have been 

1533 indicated in blue, and the corresponding 10 - 90% levels have been shown in green, (a) The 

1534 "cosine-bell" smoothing comprises 13 horizontal pixels of the streak image in this case, for a 

1535 FWHM of 0.8 1 ps, yielding a 10 - 90% rise time of 3.5 1 ps. (b) Here the smoothing comprises 

153 6 27 pixels, FWHM= 1.78 ps, and 10 - 90% rise time is 3.78 ps. "Baseline" and peak RLI are 

1537 0.185 and 1.010, respectively, in both cases. 
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Table 1 — Free Parameters in Current Model 
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Table 2 — Representative Values of Current 
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Table 3 — Optical Parameters from Streak Photographs 


Event 

Level 

Reconstructed 

RLI 



Height 

Path 

Baseline 

Amplitude 

(Fish/ 

Stroke) 


(m) 

(m) 



8715/9 

1 

30 

30 

0.0 

3 . 6 

8715/9 

2 

144 

144 

0.8 

1.5 

8715/9 

3 

274 

299 

0.5 

1.5 

8715/9 

4 , 

407 

480 

0.9 

2.2 

8715/9 

5 

488 

597 

0.6 

1.5 

8715/9 

6 

675 

864 

0.4 

1.6 

8715/9 

7 

787 

983 

0.5 

1.8 

8715/9 

8 

919 

1128 

0.5 

1.6 

8715/9 

9 

1021 

1240 

0.5 

1.5 

8715/10 

1 

30 

30 

0.1 

2.3 

8715/10 

2 

144 

144 

0.1 

1.4 

8715/10 

3 

274 

299 

0.0 

1.2 

8715/10 

4 

407 

480 

0.3 

0.8 

8715/10 

5 

488 

597 

0.1 

0.9 

8715/10 

6 

675 

864 

0.0 

0.8 

8715/10 

7 

787 

983 

0 . 1 

0 . 6 

8715/10 

8 

919 

1128 

-0.1 

0.6 

8715/10 

, 9 

1021 

1240 

0.0 

0.5 


Rise Time 

10-90% 

Error 

. (ps) 

(ps ±) 

1.3 

0.07 

0 . 9 

0.07 

1.3 

0.07 

1.7 

0.14 

1.1 


4.2 

0.21 

2.7 

0.49 

2 . 6 

0.48 
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Table 4 — Model Peak Current (kA) vs. Path Length 


Event 

Path Length (m) 

(Fish/ 

Stroke) 

0 

10 

30 

100 

300 

1000 

8715/9 

32.7 

29.3 

24.9 

19.9 

18.4 

16.8 

8715/10 

15.4 

11.7 

10.9 

10.4 

9.2 

7.8 

8725/1 

20.1 

18.4 

16.0 

12.5 

9.8 

7.1 

8725/2 

16.3 

13.3 

12.2 

10.5 

8.0 

6.3 

8725/3 

43.0 

31.3 

25.2 

23.0 

20.0 

15.8 

8725/5 

11.7 

8.6 

7.8 

7.1 

5.5 

3.4 

8726/1 

38.3 

25.1 

18.4 

17.0 

16.4 

15.8 

8726/2 

26.6 

24.8 

22.3 

18.8 

14.8 

9.9 

8726/3 

16.2 

11.8 

10.9 

10.3 

8.9 

6.2 

8726/4 

22.2 

20.7 

18.6 

15.7 

12.1 

8.0 

8728/10 

26.8 

25.5 

23.3 

18.8 

14.9 

12.0 

8728/11 

15.8 

15.2 

14.1 

11.3 

7 . 9 

5.4 

8732/1 

17.9 

17.1 

15.6 

12.5 

9.5 

7.3 

8732/2 

16.6 

15.7 

14.1 

11.3 

9.0 

6.7 

Mean Scaled 
(see text) 

1.58 

1.33 

1.17 

1.00 

0.82 

0.63 

Stdev. Scaled 

0.24 

0.11 

0.08 

0.00 

0.07 

0.12 
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Table 5 -- Comparison Among Values Averaged Over Strokes 



Path Length 

above 

Surface 

0 m 

10 m 

30 m 

100 m 

300 m 

1000 m 


10-90% Rise Time of Current Waveform: 


Average (ps) 

0.312 

0.313 

0.332 

0.460 

0.884 

2 . 18 

Std. Dev. (ps) 

0.172 

0.168 

0.158 

0.126 

0.165 

0.47 

# Samples 

24 

24 

24 

24 

24 

24 


10-90% Rise Time of Relative Light Intensity: 


Average (ps)- 


Std. Dev. (ps) 


# Samples 


Average Path (m) 


Std. Dev. Path (m) 


Effective Current-Propagation Speed: 


Average (10 m/s) 


Std. Dev. (10 s m/s) 0.30 


# Samples 24 



2.0 

3.8 

0 . 6 

1.0 

18 

20 

303 

1062 

49 

86 



1.75 

1.78 

1.81 

1.80 

0.20 

0.22 

0.30 

0.24 

24 

24 

22 

24 


2-D Average Optical-Propagation Speed: 


Average (10 8 m/s) 


Std. Dev. (10 8 m/s) 


# Samples 


Average 


Std. Dev. 


# Samples 


Peak Amplitude of Current Waveform: 




0.88 

0.79 

0.68 

0.58 

0.49 

0.10 

0 . 11 

0.08 

0.06 

0.09 

24 

24 

24 

24 

24 


Average (kA) 

22.8 

19.2 

16.7 

14.2 

11.7 

9.2 

Std. Dev. (kA) 

9.4 

7 . 1 

5.6 

4 . 6 

4.4 

4.3 

# Samples 

14 

14 

14 

14 

14 

14 


Amplitude of Relative Light Intensity: 


Average I 2 . 2 


Std. Dev. 0.9 


# Samples 14 


Average Path (m) 30 


Std. Dev. Path (m) 0 



1.1 

0.9 

0.6 

0.5 

18 

20 

303 

1062 

49 

86 
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Table 6 


Correlations with Leader Relative Brightness 


Event 

RLI Ratio 

a ( 1 0 0 ) 

ISQSsIBfl 





8715/2 


0.78 

1.06 

8715/3 


0.75 

1.20 

8715/4 


0.73 

1.23 

8715/6 


0.72 

1.08 

8715/7 


0.61 

1.40 

8715/8 


0.74 

1.16 

8715/9 

0.36 

0.61 

1.35 

8715/10 

0.10 

0.69 

1.35 

8717/3 


0.70 

1.12 

8717/5 


0.76 

1.11 

8725/1 

0.23 

0.64 

1.41 

8725/2 

0.17 

0.66 

1.32 

8725/3 

0.59 

0.54 

1.62 

8725/4 


0.76 

0.96 

8725/5 


0.64 

1 . 44 

8726/1 

0.37 

0.45 

2.63 

8726/2 

0.04 

0.72 

1.35 

8726/3 


0.65 

1.44 

8726/4 

0.00 

0.72 

1.34 

8726/5 


0.63 

1.50 

8728/10 

0.32 

0.71 

1.23 

8728/11 

0.13 

0.75 

0.98 

8732/1 

0.05 

0.71 

1.46 

8732/2 

0.08 

0 . 69 

1.34 
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Table A1 -- Camera-orientation angles, as defined in Figure A2 . 
<x represents elevation relative to the horizontal, and -y± 
represents azimuth relative to the direction to the triggering 
site. Subscript "N" stands for the NRL site, and "S" stands for 
SUNY . 


Flash : 

8715 

8717 

8725 

8726 

8728 

8732 

Parameter 







a K (deg.) 

9.98 

9.65 

10.08 

.9.96 

9.84 

9.87 

Y„ (deg.) 

1.32 

1.34 

1.97 

1.74 

1.89 

2.02 

Y s (deg.) 

-0.77 

00 

00 

T — 1 
1 

-1.61 

-1 . 74 

-1 . 16 

-1.68 
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8732/2 — Propagation Speed 
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8732/1 — Velocity Adjustment 



o 

o 

O 

O 

o 

o 

r-H 

T— H 

r-H 

r-H 

* — i 

r-H 

X 

X 

X 

X 

X 

X 

m 

in 

<N 

in 

r-H 

m 


<N 


r-H 




(s/ra) A;xoopA uopnSndojj 


Path Length (m) 












fen 



(ra/A) oSireiQ ppij 



Figure 12 
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Figure 19 
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LIGHTNING RETURN-STROKE CURRENT WAVEFORMS ALOFT, FROM 
MEASURED FIELD CHANGE, CURRENT, AND CHANNEL GEOMETRY 

J.C. Willett, D.M. Le Vine and V.P. Idone 


Return strokes (the bright channel associated with 
cloud-to-ground lightning) are the most powerful lightning 
processes. They are also responsible for most of the 
damage lightning does (to buildings, trees and aircraft) 
and for starting forest fires. Understanding the hazard 
associated with lightning reguires information about the 
current in return strokes: The information available is 

largely limited to measurements at the base of the channel 
(where it touches the ground) . Very little information is 
available about the current "aloft", which is important for 
understanding lightning process and for assessing the 
hazard to aircraft and space vehicles. 

This paper describes remote sensing research that uses the 
electromagnetic fields radiated by return strokes to infer 
the magnitude and evolution of current as it propagates 
along the channel. It consists of two parts: (a) 

reconstruction of lightning channels from stereo 
photographs; and (b) using these channels to infer the 
behavior of the current above the ground from current 
waveforms measured at the channel base and the radiated 
electric field waveforms. Among the interesting results 
are the important role channel "tortuosity" plays in the 
shape of the radiated fields and an abrupt increase in 
rise-time and decrease in amplitude just above ground. 



LIGHTNING RETURN-STROKE CURRENT WAVEFORMS ALOFT, FROM 
MEASURED FIELD CHANGE, CURRENT, AND CHANNEL GEOMETRY 

J.C. Willett 1 , D.M. Le Vine 2 , and V.P. Idone 3 

1 P.O. Box 41, Garrett Park, MD 20896 

2 Code 614.6, NASA/ GSFC, Greenbelt, MD 20771 

3 ES-351, SUNY, 1400 Washington Ave . , Albany, NY 12222 


Abstract 

Three-dimensional reconstructions of six rocket-triggered 
lightning channels are derived from stereo photographs. 
These reconstructed channels are used to infer the behavior 
of the current in return strokes above the ground from 
current waveforms measured at the channel base and 
electric-field-change waveforms measured at a range of 5.2 
km for 24 return strokes in these channels. Streak 
photographs of 14 of the same strokes are analyzed to 
determine the rise times, propagation speeds, and 
amplitudes of relative light intensity for comparison with 
the electrical inferences. Results include the following: 

1) The fine structure of the field-change waveforms that 
were radiated by subsequent return strokes can be 
explained, in large part, by channel geometry. 2) The 
average 10 - 90% rise time of the stroke current increased 
by about a factor of seven in our sample, from an observed 
0.31 ± 0.17 ]is at the surface to an inferred 2.2 ± 0.5 ps 
at 1 km path length above the surface. 3) The three- 
dimensional propagation speed of the current front averaged 
1.80 ± 0.24 X 10 s m/s over channel lengths typically greater 
than 1 km. 4) Assuming that the measured current was 
entirely due to the return stroke forced an unreasonably 
large and abrupt reduction in inferred current amplitude 
over the first few tens of meters above the surface, 
especially in cases when the leader was bright relative to 
its stroke. Therefore, a significant fraction of the 
current at the surface was probably due to the leader, at 
least in such cases. 5) Peak return-stroke currents 
decreased by approximately 37 ± 12% from 100 m to 1 km of 
path length above the surface. Because of uncertainty 
about how to partition the measured current between leader 
and return stroke, we are unable to infer the variation of 
current amplitude near the ground. 



